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FOREWORD 


This is an interim report describing the work performed during 
the period February 1982 through June 1983 under NASA Contract 
NAS8-34504, "Atomization and Mixing Study." The Rocketdyne 
project engineer for this program is Dr. Allan Ferrenberg, who 
also is responsible for the technical quality and guidance of 
the atomization work. Vance Jaqua provides technical exper- 
tise for the mixing work. Frank Kirby is the Rocketdyne Pro- 
gram Manager. Other Rocketdyne personnel supporting this pro- 
gram are Stan Pinkowski, Ed Bechtel, and Tony Exposito. This 
program was performed under the technical direction of Fred 
Braam of the NASA Marshall Space Flight Center. 


PAGE PLANK NOT FILMED 


RI/RD83-170 

iii/iv 


CONTENTS 


Summary 

Introduction 

Atomization 

Parameters Affecting Atomization 

Atomization Correlating Equations 

Droplet Measurement Techniques 

Droplet Size Measurement Results and Correlations for 

Rocket Engine Injectors 

Atomization Survey - Findings, Conclusions, and Recommendations 

Nomenclature 

References 

Mixing .... 

Injector Mixing Correlating Parameters ....... 

Miring Test Methods .... 

w — 

Cold Flow Mixing Data^ Reduction 

Review of Existing Data 

Literature Survey Conclusions and Recommendations 

Cold Flow Mixing Tests 

Liquid-Liquid Mixing Tests 

Gas-Liquid Mixing Tests 

Mixing Test Conclusions .... 

References 

Appendix A 

Appendix B 


1 

3 

5 

7 

9 

16 

22 

47 

53 

55 

63 

64 
72 
76 
80 

104 

107 

112 

125 

137 

139 

A-l 

B-l 


RI/RD83-1/0 

v/vi 


ILLUSTRATIONS 


1. Droplet Size Distributions and Relationships 

2. Droplet Cumulative Volume Distribution and Equation 

3. Derivatives of Droplet Cumulative Volume Distribution 

Data and Equation 

4. Normalized Droplet Size Distributions for Selected Injectors 

5. Effect of Gas Velocity on Droplet Size 

6. Correlation of Relative Droplet Diameter with Weber Number 

and Orifice L/d for Triplets 

7. Mixture Ratio Influence on Dropsize for Coaxial Injectors 

8. Correlation of Cold-Flow Atomization for Coaxial Injectors 

9. Influence of Combustion Gas Simulant Velocity on Mass 

Median Dropsize - Coaxial Injectors 

10. Correlation of Droplet Breakup Data coaxial Injectors 

11. Coaxial Concentric Element 

12. Impinging Triplet and Pentad Elements . 

13. Liquid/Liquid Mixing 

14. Cold Flow Gas/Liquid Mixing Measurement System 

15. Mixing Efficiency of Liquid/Liquid Triplets (ITUR 80-9) 

Depicting Elverum-Morey Factor Influence 

16. Mixing Efficiency of Liquid/Liquid Triplets (NASA CR-R-9270) 

Depicting Elverum Morey Factor Influence 

17. Mixing Efficiency of a Gas/Liquid Triplet Depicting 

Penetration Factor 

18. Normalized Cold Flow Mass Flux Profile; for LOX/GHj, 

Liquid/Gas/Liquid Triplet, 2- inch Collection Tubes 

19. Mixing Efficiency of a Liquid/Gas/Liquid Triplet Depicting 

Elverum-Morey Factor Influence .... 

20. Mixing Efficiency for Liquid/Liquid Pentads ( AFKPL- TK- 66- 147 ) 

Depicting Elverum-Morey Factor Influence 

21. Mixing Efficiency for Liquid/Liquid Pentads (NASA CR- 72827 R-8415) 

Depicting Elverum-Morey Factor Influence 

22. Mixing Efficiency for Liquid/Liquid Pentads ( AFKPL- TR- 66 1 52 ) 

Depicting Elverum-Morey Factor Influence 

23. Mixing Efficient for Gas/Liquid Pentads Depicting 

Momentum Ratio Influence * 

24. Mixing Efficiency for Gas/Liquid Pentrds Depicting 

Elverum-Morey Factor Influence 


10 

13 

14 
29 
33 

39 

44 

46 

48 

49 
67 


t O 
VO 


73 

75 


85 


86 


88 


89 

90 

92 

93 

94 

95 
9 7 


R1/RD83 170 
vi i 



25 . 


Mixing Efficiency for Gas/Liquid Pentads Depicting 

Penetration Factor Influence 

26. Mixing Efficiency for Coaxial Element Depicting 

LOX Post Recess Influence 

27. Mixing Efficiency for Coaxial Elements Depicting 

Velocity Ratio Influence 

28. Mixing Efficiency for Coaxial Elements (IL PT 73-3/) 

Depicting Velocity Ratio Influence 

29. Mixing Efficiency for Coaxial Elements (NASA CR-72/03 R-8361) 

Depicting Velocity Ratio Influence 

30. Mixing Efficiency for Coaxial Elements Depicting 

Velocity Ratio for Varying Density Ratios 

31. Mixing Efficiency for Coaxial Elements Depicting 

Falk and Burick Parameter Influence 

32. Mixing Efficiency for Coaxial Elements (IL PT 73-3 0 ) 

Depicting Falk and Burick Parameter Influence 

33. Gas/Liquid Mixing Test Apparatus 

34. Triplet #6 Mass Flux Distribution at Nominal Test conditions 

35. Triplet #6 Mass Flux Distributions at Below Nominal Mixture 

Ratio Test Conditions 

36. Triplet #8 Mass Flux Distribution at Nominal Test < onditions 

37. Liquid/Liquid Triplets - Mixing Efficiency vs 

Elverum Morey Factor 

38. Liquid/Liquid Triplets Mixing Efficiency vs 

Simulated Mixture Ratio 

39. Liquid/Liquid Triplets - Mixing Efficiency vs Momentum Ratio 

40. Liquid/Liquid Triplets - Mixing Efficiency vs 

Velocity Head Ratio 

41. Liquid/Liquid Triplets - Mixing Efficiency vs Penetration Factor 

42. Triplet No. 4 NOM MOM R 

43. Coaxial No. 5 - 43% Mixture 

44. Gas/Liquid Impinging Elements - Mixing Efficiency vs 

Simulated Mixture Ratio 

45. Gas/Liquid Impinging Elements - Mixing Efficiency vs 

Penetration Factor 

46. Gas/Liquid Impinging Elements - Mixing Efficiency vs 

Elverum Morey Factor 

47. Gas/Liquid Impinging Elements - Mixing Efficiency vs Momentum Ratio 

48. Gas/Liquid Impinging Elements - Mixing Efficiency vs 

Velocity Head Ratio 


98 

99 
100 
101 
102 
103 

105 

106 
111 
1J 5 

116 

118 

120 

121 

122 

123 

124 
128 

130 

131 

132 

133 

134 

135 


RI/RD83-170 
vi i i 


TABLES 


1. Selected Like Doublet Representative Droplet Size Correlations . 

2. Droplet Size Correlations for Triplets 

3. Droplet Size Correlations for Pentads 

4. Droplet Size Correlations for Coaxial Injectors 

5. Injector Element Correlation Parameters 

6. Summary of Previous Cold Flow Mixing Studies 

7. Representative Elements Chosen for Testing 

8. Liquid/Liquid Mixing Cold Flow Test Summary 

9. Gas/Liquid Mixing Cold Flow Summary 


24 

37 

41 

43 

66 

82 

108 

113 

126 


RI/RD83-170 

ix/x 


SUMMARY 


This report describes the results, findings, and conclusions obtained to date as 
a result of work, performed under NASA Contract NAS8-34504, "Atomization and Mix- 
ing Study," during the period February 198? through May 1983. This report con- 
tains the results of literature surveys and studies of the atomization and mixing 
characteristics of liquid rocket injectors (especially triplets, pentads, and 
coaxial injectors that may be applicable to LOX/hydrocarbon propellants), l he 
correlating parameters and other means by which such data is organized and 
reported, and the methods by which this data is obtained. Also, the results of 
the propellant mixing tests performed under this contract are reported herein. 

The major findings and conclusions of this report are summarized in the following. 

1. Tn general, very little data exists regarding the atomization and mixing 
characteristics of triplet, pentad, and coaxial elements, especially for 
impinging gas/liquid elements. Also, practically no effort has been 
expended in these areas since 1975. 

2. In general, the mixing correlating parameters provide gross estimates of 
optimum mixing efficiency, although the validity of these correlating 
parameters (and/or the validity of their reported optimum values) is 
questionable for some elements find test conditions. 

3. Perhaps the most critical parameter affecting droplet size is the local 
combustion gas velocity lield. This appears to be especially true for 
impinging elements. This is unfortunate since the actual velocity field 
in a rocket combustor and in atomization experiments is unknown. 

4. Atomization and mixing correlations and data for injector types of 
interest are presented, assessed, and summarized. 

5. The state of the art in the areas of atomization and mixing is generally 
quite poor. The physics is poorly and only qualitatively understood. 
No quantitative theories exist. The available data and correlations 
generally are of questionable validity and/or utility. Many ol the most, 
critical parameters are unknown (e.g., combustion gas velocity field. 
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multiple element effects) and/or are not simulated in tests (e.g., gas 
densities, real propellant fluid properties, combustion gas motion). In 
addition, the measurement techniques used generally employ questionable 
assumptions . 

6. The mixing characteristics of seven liquid/liquid and gas/liquid trip- 
let, pentad and coaxial, single-element injectors, representative of 
various LOX/hydrocarbon designs, have been established at various flow 
conditions . 

7. The more general objective of the mixing tests is to establish the 
existence of, and define the optimum values of, mixing correlating 
design parameters. The test results obtained to date are insufficient 
to accomplish this. Additional testing is planned. 
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INTRODUCTION 


The two major effects that reduce the performance of liquid rocket engines, 
especially larger, long-burn duration engines, are propellant mixing and evapora- 
tion. Droplet evaporation rates and chamber length requirements are highly 
dependent on the initial sizes of the droplets produced by the injectors. Hence, 
the atomization and mixing of propellants is a major concern in the design of 
rocket engines, and especially in the design of injectors. In addition, the 
potentially most powerful tools for rocket engine design are the complex computer 
codes that evaluate the combined effects of all the physical processes occurring 
in order to predict performance, stability, spacecraft contamination, heat loads, 
et. al. In order that these codes may work properly, it is necessary that the 
important physical process be quantitatively known (i.e., can be described mathe 
matically). Unfortunately, this is not the case generally, and the atomization 
and mixing processes are especially poorly known. Thus a need exists for a quan 
titative understanding, and/or extensive data, defining the atomization and mix 
ing processes as a function of injector type, propellant properties, and. opera 
ting conditions. This data consists of the mass flux distributions of the pro 
pellants and the characteristics of the droplets formed. Such information is ot 
great importance in determining critical combustion effects. Mixing efficiency 
is greatly dependent upon the initial mass flux distribution, and vaporization 
efficiency is highly dependent upon initial droplet size (and perhaps to a lesser 
extent on the initial drop velocity). Therefore, performance prediction codes 
require detailed and accurate data regarding mass flux and droplet size distribu 
tions. Also, an assessment of mixing is important in evaluating "streaking" and 
hot spots on turbine components. 

This report describes the state of the art in atomization and mixing for triplet, 
pentad, and coaxial injectors. Injectors that are applicable to LOX/hydrocarbon 
propellants and main chamber and fuel rich preburner/gas generator mixture ratios 
are of special interest. The various applicable correlating equations and param 
eters, and test data found in the literature searches are presented herein. The 
validity, utility, and important aspects of these data and correlations are dis- 
cussed. The measurement techniques employed also are presented and evaluated. 
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In addition, the propellant mixing tests performed to date under this program are 
described and summarized, results are reported, and tentative conclusions are 
presented. Additional atomization and mixing tests are planned. 

Much of the work discussed in this literature review is over 10 years old. The 
more recent atomization and mixing work has been performed as a part of injector 
or engine development programs. The primary intent of such efforts is the evalu 
ation of specific injectors operated at their design conditions, and not a deter- 
mination of the effects of the various geometric and operational variables on 
atomization or mixing. Such a "try this and see what happens" type of approach 
may be the simplest way to develop an injector, and perhaps the quality of the 
atomization and mixing data obtained with past measuring techniques did not war 
rant a detailed scientific study. However, new measuring techniques are being 
developed, and our future understanding of atomization and mixing processes 
requires a detailed, scientific study. The "try this and see what happens" 
approach is not the best, or the most cost effective, in the long run. 
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ATOMIZATION 


This section describes the state of the art in the area of liquid rocket injector 
atomization. The need for this information, the parameters of importance, the 
manners in which the data are correlated and reported, the droplet size measure 
ment techniques employed, and the specific correlations and data pertaining to 
rocket engine injectors (triplets, pentads, coaxial, and some doublet data), are 
described, discussed, and assessed. In the study of the atomization literature, 
emphasis was placed on experimental programs and empirical results directly 
related to liquid rocket injectors. The more theoretical or basic research 
efforts are to be studied in a subsequent phase of this program. A bibliography 
of the atomization reports reviewed to prepare this assessment is included in Lhe 
list of references at the end of this section. 

In liquid rocket engines, the combustion process generally is considered to be 
evaporation limited, i.e., the evaporation of the propellants is the slowest step 
in the combustion process and, therefore, very important to model correctly. 
Droplet evaporation rate is a strong function of droplet size and velocity rela 
tive to the gas phase. Some computer codes calculate drop velocity and motion. 
This is important in properly assessing evaporation, stability, spacecraft con- 
tamination by ejected propellant droplets, performance, and wall effects (e.g., 
wall film buildup, heat transfer). Droplet acceleration is due to an imbalance 
between droplet inertial forces (a function of drop diameter cubed) and drag 
forces (generally a function of drop diameter squared). Thus, the droplet size 
is an important parameter in the assessment of droplet evaporation rate and 
motion. 

It is possible to write the equations governing the motion and evaporation of a 
droplet. The forms of these equations and most of the parameters are known 
fairly well over many operating regimes of interest. The equations are ordinary 
differential equations that can be solved readily. However, as with all differ 
ential equations, any such solution requires knowledge of the initial or boundary 
conditions. And this is the problem these initial conditions are not known well 
enough. These conditions are the droplet size and velocity distributions at the 
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locations where the droplets are formed. Given these initial conditions, the 
governing equations can be solved, and this is precisely what the combustor codes 
do. However, errors in the initial conditions produce corresponding errors in 
the predictions. 

This problem has long been recognized and a number of experimental programs have 
been performed to establish these initial conditions. Due to the complexity of 
the physical processes occurring during atomization, these initial droplet condi- 
tions generally are characterized by empirical correlations. Some of these cor 
relations and experiments are described later. Both mixing and atomization 
experiments often are performed with propellant simulants. This introduces a set 
of corrections that must be employed to extrapolate to the actual propellants of 
interest. Another set of corrections generally must be applied to extrapolate 
the test data to the operating conditions (pressures, temperatures, etc.) of 
interest. Thus, the establishment of these critical initial conditions depends 
entirely upon a relatively small quantity of empirical data, relating the effects 
of a few of the many parameters affecting these complex physical processes; and 
several sets of corrections to this data. 

The utilization of such atomization data and correlations in the combustor analy 
sis codes is a major source of difficulty and error. This has been demonstrated 
repeatedly in code development programs at Rocketdyne. The three major perform- 
ance codes in use at Rocketdyne (TPP, CICM, and SDER) all attempt to use such 
correlations. In all three cases, it has been found to be necessary to modify 
the experimental correlations to force the codes into agreement with large-scale, 
rocket engine performance tests. Such "calibration" of computer models with the 
actual hardware they attempt to model is a standard procedure when dealing with 
complex unknown phenomena, although it is obviously a technique of last resort. 
Codes that are calibrated in such a manner can be relied upon to produce good 
results as long as they are applied to designs and conditions not significantly 
different from those that they were calibrated against. However, the accuracy of 
the codes becomes increasingly questionable as they are applied to situations and 
problems significantly different from the calibration points. 
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The poor performance of these drop size correlations indicates that something is 
wrong with the measurement techniques, the correlations developed from the mea 
surements, and/or the manner in which they are applied. The assessment of the 
state of the art, as described in the remainder of this section, provides reason 
to suspect all of these. 

PARAMETERS AFFECTING ATOMIZATION 

The physical processes occurring during atomization cannot be reduced currently 
to sets of equations derived from basic physical principles. The most common 
case of the break up of a single jet of liquid has been theoretically studied for 
over 100 years, and these theoretical studies have been unable to predict, to an 
adequate degree, the characteristics of the droplets produced. Impinging streams 
and other fan-forming injectors also have been investigated theoretically. These 
studies and experimental efforts, combined with the strictly empirical invest iga 
tions of others, provide an indication of the parameters of importance in the 
atomization process. However, there is considerable disagreement regarding the 
relative importance of specific variables. 

The properties of a liquid propellant that are considered of importance are the 
surface tension, viscosity, and density. For a propellant injected as a gas, the 
only thermodynamic property generally considered of importance is the gas den- 
sity. The geometric variables of importance for impinging- type injectors are the 
orifice diameter, oiifice length, orifice entrance conditions, number of orifices 
(triplet or pentad), free jet distance (i.e., distance from the orifice to the 
impingement point), and impingement angle. Flow variables to be considered are 
the velocities of the liquid streams and the existence of turbulence in these 
streams . 

For coaxial injectors, the geometric variables of potential importance are the 
propellant flow areas, the inner tube (LOX post) wall thickness, and the recess 
of the LOX post. The flow variables of greatest concern are the liquid velocity 
and the relative gas to liquid velocity. 
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Another parameter that has been shown to be of very great importance is the 

velocity of the combustion gases relative to the injected fluids. This parameter 
affects the aerodynamic breakup of the large droplets and ligaments after they 
have separated from the spray fan. This is often referred to as secondary atom- 
ization and many basic research efforts have been performed to evaluate this. 

Separation of atomization into primary and secondary atomization processes is 
certainly an oversimplification, but it has been employed. The importance of the 
combustion gas velocity is unfortunate, since the actual velocity field in the 
combustor cannot be determined adequately. The combustion gas velocity field 
depends entirely upon the droplet evaporation rate and distribution, which in 

turn is highly dependent on initial droplet size, which in its turn is greatly 
affected by the combustion gas velocity field. Thus, all of these phenomena are 
interrelated highly and must be considered together. Even in cold flow tests, 
the local gas velocity field in unknown. There is no such thing as "spraying 
into still air," as the spray itself transfers momentum to the gas and sets it in 
motion . 

Another important parameter that is difficult to quantify is the liquid veloc 

ity . Generally, this is assumed to be the average velocity at the orifice exit 
assuming the orifice is flowing full. But this is not the velocity of the liquid 
in the fan, or of the ligaments, which is probably the velocity of greater 
concern. 

The effects of combustion on atomization are unknown. Matching the density and 
velocity field of the combustion gases in a cold flow atomization test may not be 
sufficient. Burning droplets may break up differently (secondary atomization) 
than nonburning droplets due to the effect of the burning gas envelope about them. 

Our knolwedge of many of the parameters affecting rocket engine atomization comes 
primarily from the study of doublets. Details of many of these studies are 
included later in this report where the relative importance of the parameters 
affecting atomization are discussed in greater detail. 
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ATOMIZATION CORRELATING EQUATIONS 


The objective of atomization studies is to develop quantitative relationships 
defining the effects of the various governing parameters on the characteristics 
of the spray. The characteristics generally most desired are a representative 
droplet size (a mean or average value) characterizing the spray, and a droplet 
distribution defining how the number of droplets in the spray varies with droplet 
size. Other parameters of interest are the distribution of droplets in space, 
and the velocities of the droplets. 

The data generally most desired is the number of drops of each size in a given 

spray. Often the data is obtained in the form of numbers of droplets counted, n, 

in each of many uniform size ranges (e.g., 5 to 10u, 15 to 20yi, etc.) as 

shown in Fig. la. However, this discrete form of data presentation has the 

undesirable characteristic that as the width, AD, of the size ranges is varied 

the count will change. In order to quantify the data in the form of a continuous 

mathematical expression, the data is often converted to the form of Fig. lb. 

Here, the number of droplets per unit size range (e.g., per micron) is plotted. 

This continuous function is called a distribution function f, and is determined 

by evaluating n ~ ^ as AD approaches zero. n(D) is the number of drop 

AD * « AD 

lets having diameters between D - “ and D + This distribut ion 

function is the mathematical expression that best defines the size distribution 

of the droplets . All the other forms and techniques for expressing droplet size 

distributions can be derived mathematically from this distribution function, f, 

where 

f - 

AD-* 0 

Another useful function is the fraction o( droplets in the spray at diameter D, 
which is 


f/n t 


1 1m 
AD-»0 


n/n t 

AD 


where n^. - total droplet count 
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Figure 1. Droplet Size Distributions and Relationships 
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In some applications, the volume of liquid in the spray as a function of drop 
diameter is of interest. Multiplying the number of drops in each diameter range 
by the volume of a drop of that diameter converts the number distribution to a 
volume distribution (Fig. lc). It is also possible to construct a volume distri 
bution function (not shown in Fig. 1) analagous to the size distribution func 
tion. Another useful representation is a cumulative distribution. The cumula 
tive number, N, of droplets at any diameter D is the sum of the number of all 
droplets having diameters less than D (Fig. Id) . 

0 D jij 

N(D) = 2n = [ f dD or (t») - f (0) 

ok 0 

Similarly, the cumulative volume distribution, V, of the droplets at any diameter 
D is the sum of the volumes of all drops in the spray having diameters less than 
D (Fig. le). 


0 

V( 0) = £ 
0 


v - 


Jo 


11m 

ftD-»o 




dD 


or 


'dD' 


D 


1 

6 


D 3 f 


The normalized cumulative volume distribution, R, (Fig. If) is the volume distri 
bution divided by the total volume of all the droplets measured, i.e.. 



The cumulative volume distribution (often normalized) is the most commonly util 
i zed manner for graphically presenting the data. 

The mathematical expressions defining the drop size distribution that are 
encountered most commonly are 

= AD 5 exp ( -BD n ) Nukiyama - Tanasawa 

jjr = BnO* 1 " 1 exp ( BD n ) Rosin Rammler 

au 

~ = * ff“ 1/2 exp (-6 2 Y 2 ) Log probability 

where Y = In (D/D) 
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and A, B, n, and 6 are the adjustable constants. Many other distribution 
functions have been utilized and a more complete list and description of these is 
contained in Ref. 11. Since: 


dV 1 n 3 
do * 6 ’ 0 


dN 1 „3 

dO ' 6 * ° 


f 


these distributions actually define the desired distribution function, f. Given 
the distribution function, f, the cumulative volume (or number) distribution can 
be obtained by integrating the distribution function, f, over various size ranges. 


More often, the data is plotted in terms of the cumulative volume, or normalized 
cumulative volume, versus drop diameter. For some droplet measurement tech- 
niques, particularly the frozen wax technique, it is this cumulative volume (or 
mass) distribution which is measured directly. Evaluation of the slope of this 
distribution then can be performed to define the droplet number distribution or 
distribution function. Cumulative distributions tend to "smooth’' the data, mask, 
inaccuracies due to too few droplet measurements, and reduce the apparent differ- 
ences between different distributions. Usually, only the cumulative volume dis 
tribution is reported, so this problem is overlooked often. One very comprehen- 
sively reported investigation (Ref. 66), which presented all of the data and 

plots of V and ^ versus D, demonstrates this problem. Figure 2 is a cumu 
r dD 

lative volume plot of the data from one set of droplet size measurements as pre 
sented in that report. The data appear to be in good agreement with the integral 
of the particular distribution function chosen. However, the plot of the actual 

and the data (Fig. 3) demonstrate that this appar- 
distribution function and the data is misleading. 


volume distribution, — , 
ent agreement between the 


I 
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Figure 2. Droplet Cumulative Volume Distribution and Equation (Ref. 66) 
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Figure 3. Derivatives of Droplet Cumulative Volume Distribution 
Data and Equation (Ref. 66) 
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Often the drop distribution is characterized by a single value- a representative 
drop diameter. Some such common representative diameters are the mean or average 
drop diameter, volume mean diameter, Sauter mean diameter, or mass median diam 
eter. The mean or average diameters are defined according to the following 
general relationship: 


0 


pq 


1 




p-q 


2 n . 


D 


q 

i 


where i denotes size range considered. 


n = number of droplets in size range i 
i 

D = middle diameter of size range i 
i 

Thus, D 1q is the linear average diameter of all the droplets measured, is 

the diameter of the droplet having the average volume of all the droplets meas 
ured, and D 32 is the diameter of the droplet whose volume to surface area ratio 
is the average of all the droplets in the spray (referred to as the Sauter mean 
diameter). The mass median droplet diameter is the droplet whose size is such 
that one half of the mass (or volume) of the spray is contained in droplets 
having a larger diameter. On a plot of K versus D (Fig. If), the mass median 
diameter is the diameter occurring at a value of R ~ 0.5. All of these repre 
sentative diameters can be calculated from t he distribution function, f. 

Most of the correlations developed to define the effect of various geometric and 
operational variables on droplet size define this effect in terms of the intlu 
ence these parameters have on one of these representative diameters. These 
usually take a form: 

Representative diameter = A ,n ^ X-j 
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where the X terms are the variables, or collections of variables, of interest, 

and A, m, p, are the adjustable constants by which the relationship is 

made to fit the data. In some instances, sums of terms, each similar in form to 
the right-hand side of the above equations, are employed. It must be recognized 
that such relationships as this do not completely characterize the spray, and 
that two sprays with the same mass median or Sauter mean diameter are not the 
same. It is often the smallest drop size and/or the largest drop size that are 
of greatest importance (e.g., in assessing stability and performance respec- 
tively). The mean or median droplet size does not provide this information. 
Thus, it is important to also characterize the droplet size distribution- 1 .e. , 
to establish the correlating equation defining the distribution function, f. 


DROPLET MEASUREMENT TECHNIQUES 

A variety of techniques have been employed to measure droplet sizes. All of 
these techniques are subject to inaccuracies and questions associated with the 
basic assumptions employed, their manner of use, and/or the quantities of data 
usually obtained. Details of these techniques can be found in the literature, 
including some of the references contained in the bibliography. The discussion 
here is limited to the three primary techniques previously employed to obtain 
atomization data for rocket engine injectors. The findings obtained regarding 
rocket engine injection atomization utilizing these techniques is discussed 
subsequently . 


Imaging Techniques 


These include photography and holography and have been the most extensively 
employed methods for droplet sizing. They generally require a fairly dilute 
spray and offer the advantage of actually "seeing" the droplets as they exist at 
the point and time where knowledge of their size is desired. Although multiple 
exposure techniques can be employed to obtain droplet velocity data, none of the 
experimental programs discussed herein did so. As will be discussed shortly, 
velocity information is essential to the determination of accurate droplet size 
distributions when imaging techniques are employed. Imaging techniques have been 
employed to measure droplet sizes in reacting flows. This is an important and 
valuable feature that apparently has only been employed for the case of a rocket 
engine combustor in the investigations reported in Ref. 26, 72, 73, and 76. 
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A major problem with the use of imaging, techniques has been the need for human 
analysis of the images. Although computerized techniques have been developed 
recently for analysis of photographic images, all of the rocket injector atomize 
tion work described herein employed at least some degree of human involvement, in 
the analysis of the droplet images. It is necessary for someone to determine 
which droplets are to be counted (i.e., which droplets are in focus), and in most 
cases, to manually measure the droplet sizes. This causes errors in two ways 
human judgment and insufficient droplet counts to define the spray. 

Another problem associated with imaging techniques is the time (i.e., cost) 
required to manually identify, count, and measure the droplets. This often pre 
vents the counting of a sufficient number of droplets to assure an accurate dis 
tribution. A large number of droplets must be counted. The number of small 
drops may be over 1000 times as great as the number of large ones, and yet these 
large drops are often the most important to include. In Ref. 5, it is calculated 
that it is necessary to count 5500 droplets to be 95% confident that the Sauter 
mean diameter is correct to within 5%. Rarely are so many droplets counted per 
sample with imaging techniques. 

Perhaps the most important problem associated with imaging techniques is that 
these techniques only measure the concentration of droplets in a given volume of 
space (i.e., spatial distribution) rather than the true droplet distribution, the 
temporal distribution. This problem is recognized in the older literature 
(Ref. 11), but appears to have been neglected by many others. The nature of this 
problem is discussed in detail in the following. 

In a steady-state flow of droplets, the number of droplets, and the number of 
droplets of each size entering a particulat region in space per unit time must, be 
constant. It is possible to write a droplet number conservation equation (anal- 
agous to a mass conservation equation) as follows: 

N 1 = '1 A V 1 

where ^ ^ number of droplets of size group 1 entering a region or control 
volume (drops/sec) 

p. = local concentration of droplets of size 1 (drops/cm ) 
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A = cross sectional area of region perpendicular to the direction of 

2 

flow (cm ) 

= velocity of the droplets of size group 1 (assuming all drops In 
size group 1 are travelling at the same velocity) (cm/sec) 

Now, the temporal distribution of droplets produced by an injector can be 

obtained by counting the droplets per second of each size group i crossing A, 

that is, by measuring the values of the terms. As long as the droplets are 

moving in one direction (i.e., the spray is not spreading), measurements of the 

N values at any location in the spray will not change. However, the imaging 
i 

techniques measure the droplet concentrations, i.e., the p. terms. As long 

as the droplet velocities remain constant as the spray moves downstream, these 

p terms also will remain constant. However, if the droplet velocities 

i 

change in such a way that the smaller drops are no longer moving at the same 
speed as the larger drops, then the Pj terms also change. The values 

must remain constant for this is a steady flow situation. Thus, an imaging tech 
nique measures the p. terms, and the ratios of the p^ terms is not the 
ratios of the actual number of droplets of each size group in the spray. The 
only time that the imaging techniques produce true drop size distributions is 
when all the droplets move at the same velocity. This condition rarely, if ever, 
occurs in nature or in experiments. 

One particularly noteworthy effort that appears to demonstrate this effect is the 
work of George (Ref. 72. 73. 76). Tn these experiments, measurements were made 

at several axial locations downstream of the injector utilizing an imaging tech- 
nique (holography). In all these tests, the gas velocity exceeded the liquid 
injection velocity. In such a case, the small droplets would be accelerated more 
rapidly than the larger droplets. This would cause the spatial concentrations of 
smaller drops to decrease faster than for the larger drops as we move downstream 
from the injector face. Thus, we should expect to see more larger drops than 
smaller drops in the holograms as we move downstream. This effect was observed 
(Ref. 73) and was quite significant. The value of D was found to increase by 
50 microns or more over a 2 inch change in axial distance. Also, a simple com- 
puter simulation of droplet dynamics in a constant velocity gas flow (Ref. 84) 
demonstrates significant differences (40% or more variation in representative 
droplet sizes) between temporal and spatial measurements. 
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This problem of spatial versus temporal distributions places some doubt on the 
utility of the results obtained with imaging techniques. Not only are the dis- 
tributions measured not the true distributions, but the distributions will vary 
at different locations due to differences in velocities of the various size drop 
lets. This may account for much of the disagreement between various investiga- 
tors. The only situation in which the temporal and spatial distributions are the 
same is when the velocities of the droplets are not a function of droplet size. 
Spraying into "still" air will never produce this condition; and spraying into 
flowing air only will approximate this condition beyond some unknown distance 
downstream from the injector (where the droplets and gas velocities are equal).lt 
is the temporal distribution that is needed for the combustor models. 

Liquid Droplet Capture Technique 

This technique involves the capture of a sample of the spray on a solid surface 
(e.g., a glass slide) or in another liquid. The droplets captured are measured 
under a microscope or photographed for later analysis. Most of the work utiliz- 
ing this type of technique was performed before 1960 and the technique seems to 
have been supplanted, to a large extent, l>y photographic and other methods. In 
many cases the captured droplets are no loi ger spherical (e.g., droplets captured 
on a surface) and corrections to account fox’ this must be applied. 

This type of measuring technique requires the use of highly nonvolatile liquids 
and, when the droplets are captured in another liquid (e.g., a heavy oil or gly 
cerine) , further requires that the droplet liquid be immiscible in the capturing 
liquid. This limits the choice of liquids that can be utilized. Also, the drop- 
lets must be captured gently so as to prevent droplet shattering. 

In many applications of this technique, it is obviously the temporal distribution 
of droplets that is obtained. For some sampling methods, however, there is some 
question as to whether it is the spatial or temporal distributions that are 
measured. Such questionable methods include the slide "waving" technique, where 
a glass slide is passed rapidly through a spray. 
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Like the imaging technique, the liquid droplet capture technique requires consid 
erable manpower to count and size the droplets. Thus, the size of the sample 
counted may be a serious source of error. Also, this technique requires the 
spray to be diluted sufficiently to prevent a significant amount of coalescense 
of the droplets in the sample. In order to accommodate this requirement, one 
technique often employed is a spray splitter. The spray impinges on a plate con- 
taining a hole or slit through which only a portion of the spray may pass. Only 
this small portion is allowed to fall on the collection surface. This same pro 
cedure also is used occassionally with imaging techniques to dilute the spray. 
One aspect of this spray splitting procedure that occassionally is overlooked is 
the effect and probability of droplets colliding with the edge of the splitter 
plate. Such collisions can shatter droplets thereby causing the sampled spray to 
have a droplet distribution different from that of the main spray. This problem 
is analyzed in some detail by Dickerson (Ref. 47). 

Droplet Freezing Techniqu e 

This technique has been applied extensively in the study of rocket engine injec 
tors. Much of this work that is related directly to rocket engine injectors was 
performed at Rocketdyne during the period 1967 through 1975, and utilized wax as 
the injected liquid. Fluids other than wax have been used and droplet capture 
and freezing in liquid nitrogen also has been performed. All of the work 
described herein utilizing this technique was done with wax. 

The frozen wax technique offers several advantages over other methods. The 
liquid wax is injected into the atmosphere or a large pressure vessel where the 
droplets rapidly cool and solidify, and then are collected and sampled. The sam- 
ple then is subjected to a sieving operation where the wax droplets are separated 
into size groups. Each size group then is weighed and a plot of droplet mass 
(volume) versus size is constucted. Thus, the cumulative volume, volume distri- 
bution, and mass median diameter are measured directly, without the great time 
and manpower associated with the sizing and counting of individual droplets. 
Also, the true or temporal distribution of droplets is measured, since all the 
droplets produced by the spray over a long period of time (several seconds) are 
collected. And finally, the number of frozen wax droplets included in the sample 
is on the order of millions. This technique does not suffer from a lack of a 
sufficient sample size to be statistically accurate. 
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One serious disadvantageous feature of the hot-wax technique is the limited 
choice of materials that can be applied conveniently. Since the properties of 
the actual propellants are different from the simulants, it is necessary to 
establish the effects of these properties (surface tension, density, and viscos 
ity) on the atomization process. Thus, the capability to perform tests with dif 
ferent fluids having widely varying properties is important. Another feature of 
the wax technique that merits consideration is the density increase upon freez- 
ing. Because of this, some earlier investigators have corrected the wax droplet 
sizes by multiplying the measured droplets’ volumes by the ratio of the solid to 
liquid densities. However, the physics of the freezing phenomena indicates that 
the droplets should freeze on the outside first. If this is correct, then the 
frozen drops should be hollow and no density change drop size correction is 
required. Dickerson (Ref. 47) has discovered that the droplets indeed are 
hollow, and that the volume of the central void is approximately equal to the 
size change due to freezing- at least for the larger drops. 

One of the most serious criticisms of most of the hot-wax investigations involves 
the problem of defining the temperature (and hence the properties) of the liquid 
wax during atomization. In most investigations, the hot liquid wax is injected 
into a relatively cold gas (e.g., the atmosphere). For these cases, it is neces 
sary to question whether the wax has cooled significantly prior to the completion 
of atomization. Zajac (Ref. i>8) presents data showing that the surface tension 
and viscosity of the particular wax utilized (shell 270) increase by 12% and 83%, 
respectively, between 93 C (the nominal injection temperature) and 66 C (slightly 
above the wax fusion temperature). Certainly, the surfaces of the wax ligaments 
and droplets must be cooler than the bulk wax injection temperature. Thus, the 
wax properties at the injection temperature may not be the same as those existing 
during atomization. Longwell (Ref. 1) presents results suggesting the wax tech 
nique erroneously may give large droplets due to viscosity increases as the 
liquid cools during atomization. However, Hasson and Mizrahi (Ref. 23) present 
extensive data demonstrating that the wax technique produces significantly and 
erroneously small droplets (they corrected for an assumed shrinkage of the drop 
lets upon freezing, but this correction is not great enough to account for the 
observed difference). Several investigators (Ref. 29, 70, 71, 78) performed hot 
wax experiments in which the wax was injected into hot gas and subsequently 
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cooled after atomization was complete. In these investigations the potential 
problem of wax cooling during atomization should have been eliminated or 
minimized . 

DROPLET SIZE MEASUREMENT RESULTS AND CORRELATIONS 
FOR ROCKET ENGINE INJECTORS 

This section presents all the pertinent atomization results that were found 
relating to triplet, pentad, and coaxial injectors. Since very little data per- 
taining to these injector types has been found, and since most of our knowledge 
of atomization comes from the study of like doublets, a discussion of like doub 
lets also is included. 

The expressions relating representative droplet diameter to injector geometry, 

operating conditions, and environment vary with each investigator. The most 
common representative diameters utilized are the Sauter mean ^^ 32 ^ ’ volume or 
mass mean (D 3Q ), and mass median (5). Conversion between these diameters 
requires that the droplet size distribution be known, and a generalized conver- 
sion requires that the distribution function be known and integrable. Generally, 
such information is not available, so a direct comparison between these repre 
sentative diameter equations cannot be accomplished (one exception to this is 

described later). However, inspection of the exponents of some of the more 

important variables (e.g., liquid velocity, and orifice diameter, d j ) , 

indicates considerable disagreement between these equations. This may be due to 
the previously discussed questions and problems regarding the measurement tech 
niques, testing over different ranges and conditions, the use of different 

fluids, unmeasured and/or uncontrolled variables (the most important being the 
local gas velocity, V ), and/or other unknown causes. In a few cases, these 
drop size equations contain variables that are not varied significantly during 
the testing. In many cases (all cases for the triplet, coaxial, and pentad 
injectors), all of the potentially significant variables have not been tested. 
The equations developed from such data are incomplete. All of the droplet size 
equations described herein are strictly empirical or are based only in part ot. 
very limited theoretical considerations. 
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The most extensively studied rocket engine type injector is the like doublet. 
Much can be learned through a review of these doublet studies about atomization 
processes and measurements in general that is relevant and important to the study 
of triplets, pentads, and even coaxial injectors. Table 1 is a list of selected 
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Furthermore, in an earlier investigation (Ref. 13), Ingebo established a rela 
tionship for the effect of injected liquid properties on droplet size as: 


.25 

V 

°30 “ V 

However, this relationship was determined in experiments involving crosscurrent 
injection of single streams into flowing gases. Its applicability to cocurrent 
injection may be questionable, yet, it is often utilized. 


In a subsequent study (Ref. 26), Lngebo utilized a moving camera to photograph 
burning ethanol droplets and measured their velocities. At a distance of 0.1 m 
from the injector, the drops were observed to be traveling at a higher velocity 
than their injection velocity. Most importantly, the small droplets were 
observed to have been accelerated much more than the larger droplets. 35 micron 
droplets had undergone a velocity increase 9 times as great as 344 micron drop- 
lets. Again, this would indicate that the spatial concentration of each size of 
droplets would be rapidly and differently varying with distance from the injec 
tor. The effect this would have on the measured spatial droplet size distribu 
tion apparently was not considered. 


In 1964, several studies (Ref. 36 through 38) were reported by investigators at 
Aerojet General. Brown (Ref. 38) captured droplets on glass slides that were 
produced by the injection of a stream of liquid into flowing cold and hot (up to 
nearly 1300 K) gas. One of the important features of this work was the recogni- 
tion that the spray affects the gas velocity. An attempt was made to quantify 
this effect in a droplet size relationship with a term containing the mass flow- 
rate ratio of liquid to gas. In another of these investigations, Wolfe and 
Anderson (Ref. 37) performed experiments and developed a relationship for the 
breakup of large droplets (i.e., secondary atomization) based upon the earlier 
work of Weiss and Worsham (Ref. 29). This relationship includes a liquid proper- 
ties effect of the form, 
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In the early 1970s, photography and the new technique of holography were utilized 
at AFRPL. to measure droplet sizes (Ref. 54, 66, 72, 73, 76). Kuykendal (Ref. 54) 
investigated the effects of liquid velocity, orifice diameter, impingement angle, 
stream alignment, orifice length, and surface finish for like doublets flowing 
water. Droplet size equations were developed to define these effects, but the 
average drop sizes were based upon a relatively small drop count (occasionally 
less than 100), and these equations appear to disagree greatly with most other 
similar studies. George (Ref. 72, 73, 76) utilized holography to measure droplet 
sizes in both hot flow (hydrazine drops burning in nitrogen tetroxide) and cold 
flow (water in N 2 ). The form of the droplet size correlations developed in 
that effort are presented in Table 1. 

During the late 1960s and early 1970s, a very elaborate hot-wax capability was 
developed at Rocketdyne and many atomization investigations were performed (Ref. 
47, 48, 50, 51, 52, 55. 56, 58, 60, 61, 63, 64, 65, 67-71, 74, 78). Two of the 
most comprehensive of these investigations are those of Dickerson, et . al . (Ref. 
47) and Zajac (Ref. 58). Their correlations for like doublets are presented in 
Table 1. Both of these efforts were performed by spraying wax into "still" air. 
Dickerson's droplet size correlations, as reported in Ref. 47, are not in 
agreement with his subsequent paper (Ref. 52). Discussions with Dickerson 
revealed that the earlier liquid velocity data was incorrect, and the correla- 
tions of his latter paper include the correction. Dickerson evaluated the atom 
ization characteristics of a variety of impinging injectors, with great emphasis 
on doublets. Experiments were performed with impinging fans from unlike pairs of 
like doublets utilizing water as the other fluid. These tests indicate that 
impinging fans tend to broaden the droplet size distribution but have little 
effect on D. Droplet size distributions for several of the injectors tested are 
presented in Fig. 4. Note that both axes have been normalized in such a way that 
all distributions must pass through the point (1.0, 0.5). Also, as previously 
discussed, it is the slope of these cumulative volume distributions that truly 
defines the spray. Thus, the apparently small differences in the plots of 
Fig. 4. are, in fact, large differences of great importance. 
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The single most comprehensive study of atomization of rocket engine-type injec 
tors is the work of Zajac (Ref. 58). Zajac examined the effects of liquid veloc- 
ity, orifice diameter, velocity and diameter ratios, orifice length, free jet 
length (distance from orifice to impingement point), angle of impingement, ori- 
fice entrance conditions (geometric and flow conditions), misimpingement , and 
propellant miscibility for like and unlike doublets are well as triplets and pen- 
tads. In addition, he measured transient pressure distributions within the free 
streams (a measure of velocity profile and turbulence). Zajac found that streams 
flowing turbulent acted considerably different than laminar streams with regard 
to atomization (Dickerson had neglected this, but earlier investigators, e.g., 
Dombrowski (Ref. AO), already had indicated this). Thus, it was necessary to 
establish two droplet representative diameter equations, one for turbulent and 
one for laminar. Velocity profile also was found to be important, but only in 
laminar flow. Free stream breakup prior to impingement was shown to be important 
and can occur at a free stream length of from 5 to 10 orifice diameters in turbu 
lent streams. The much higher gas densities in a real combustor could cause 

breakup in shorter lengths. 

The state of the art circa 1971 was that the wax technique yielded sufficient 
quantity and apparent quality of data to define droplet sizes and size distribu- 
tions of hot wax droplets sprayed from like doubltes into still air. Several 
problems remained, as follows: 

1. How valid is the hot wax technique? Does wax significantly change prop- 
erties before atomization is complete? 

2. How can hot-wax results be correlated to that of real propellant? 

3. What is the effect of the actual rocket combustor environment (hot, 
high-density combustion gases moving at high velocity) on the atomiza- 
tion process? 

In an attempt to solve some of these problems, tests were performed with combine 
tions of waxes to examine viscosity effects (Ref. 69) and a large pressure tank 
was utilized to simulate high-density gases (Ref. 64, 69, 74, et. al.). In addi- 
tion, several attempts (Ref. 65, 68, 74, et. al . ) were made to validate these 
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droplet size correlations by utilizing them in computer models of rocket engines 
and comparing the results of these models with the actual hot-fire tests the 
models were attempting to simulate. In one program (Ref. 65), a test engine 
operating on wax and liquid oxygen was utilized. Although all of these efforts 
reported some degree of success, these three basic questions still remain essen 
tially unanswered. 

One of the most unknown aspects of this problem was (and is) the effect of the 
combustion gas velocity on droplet size. The actual velocity field existing in 
and around the spray in cold-flow experiments is never measured. The actual vel- 
ocity field existing in and around the spray in an operating engine also is 
unknown. And finally, the effect of a known flowfield on the formation and 
breakup of a spray fan or stream (primary atomization) is essentially unknown. 
There is, however, a considerable body of work performed to evaluate the effects 
of gas flowfields on the deformation and breakup of individual droplets (second- 
ary atomization) . Such efforts demonstrate the great complexity of this latter 
process. 

In an effort to establish the effect of gar velocity on the size of droplets pro 
duced by impinging liquid streams, experiments have been performed in low 
pressure wind tunnels. In such experiments, the gas velocity is defined as the 
velocity that existed prior to the introduction of the spray. The effect of the 
spray on the gas velocity, although often recognized, is not taken into account- 
very crudely included by Brown (Ref. 38), and is not measured. Similarly, the 
liquid velocity in the gas is assumed to be the average liquid velocity at the 
injector orifice exit, and not the actual liquid velocity in the spray fan. 
Thus, in attempting to correlate this very important effect of relative gas vel 
ocity (gas velocity relative to liquid velocity), the velocities used are incor- 
rect and are, at best, only representative of the true velocities. Despite this, 
these experiments do provide an indication of the importance of the relative gas 
velocity. Probably the most extensive of these efforts for impinging like doub- 
lets are the work of Ingebo (Ref. 15) and George (Ref. 76), as previously dis- 
cussed, and the latter investigations of Zajac (Ref. 70 and 71). 
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Zajac utilized a doublet injecting hot wax cocurrently into a ducted hot 
60 C) nitrogen stream. He separated the atomization process into two parts 
(i.e., primary and secondary) and studied these separately. In the primary atom- 
ization study, the effect of constant velocity and accelerating gas streams on 
the sizes of droplets initially formed was investigated. In the secondary atom- 
ization study, known droplet size distribution sprays were subjected to accelera- 
tions to observe droplet breaking. The rate and degree of acceleration was con- 
trolled by varying the length and area of the duct downstream of the injection 

location . 

Zajac found that many of the parameters investigated in his previous work were of 
little importance compared to the effect of relative gas velocity. Much of his 
data was plotted in the form of Fig. 5, showing droplet size versus a nondimen 
sionalized relative gas velocity. Note that all of the investigations in which 
the liquid was injected into "still" air would be plotted at the -1 value of the 
nondimensionalized velocity. Shown in Fig. 5 are volume mean diameter data from 
Ingebo showing the effect of gas velocity on droplets produced by two different 
injectors, mass median diameter data from Zajac, and the calculated droplet size 
based upon tests with v g = °- The data from Za J ac Presented here was obtained 

with a constant gas velocity (i.e., duct area remained constant). Figure 5 
demonstrates the great efect of gas velocity on droplet size. 


Based upon his experiments with accelerated and constant velocity flows, Zajac 
constructed the droplet size correlation equations shown in Table 1. These 
equations compute the mass median droplet size based upon the gas and liquid vel- 
ocities, the droplet size, "5, occurring when the maximum gas velocity equals 
the injected liquid velocity (i.e., V g -V L /V L = 0), and a parameter , which 

includes the distance over which the gas is accelerated (at V & = constant, L 
infinity). The parameter D c is computed from the liquid velocity and orifice 

diameter. A study of the derivation of these equations indicates they are 
applicable to turbulent flow only. 


The correlation of Zajac is, to some extent, supported by the earlier work ol 
Ingebo, and the very important effect of relative gas velocity is demonstrated. 
Unfortunately, the application of such results to real combustors is difficult 
since the combustion gas velocity cannot be defined adequately. 
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Figure 5. I'ft'eet of Gas Velocity on Droplet Size 





Since about 1975, there has been very little atomization work directly relateable 
to rocket engine like doublet injectors. This is certainly not because the prob- 
lem has been considered solved. Despite all the earlier efforts to define the 
initial droplet sizes produced by like doublets in combustors, our knowledge in 

this area is very crude and/or qualitative. All of the droplet size data to date 

is of questionable accuracy and/or validity due to real or possible droplet size 

measurement technique problems as previously discussed. The droplet size correl- 

ations and distributions developed from this data are generally, strictly empiri- 
cal. They are mere curve fits of the test data and, as such, may be neglecting 
important untested variables and are certainly not of the proper form. These 
correlations are based upon data that demonstrated poor or usually unknown 
repeatability, considerable spread, and often a relatively low quantity of drop 
let counts. To some extent, these features of the data are masked by the exten- 
sive use of semilogarithmic plots of the data and cumulative droplet size distri- 
bution plots. 

Perhaps the greatest problems involve the application or utility of the atomiza 
tion data. Extrapolation of the cold- flow data using wax or other liquids to the 
actual propellants and to the conditions existing in a rocket combustor requires 
many questionable assumptions and estimates. One of the most important and, 
unfortunately, most questionable of these extrapolations involves the combustion 
gas velocity, as previously discussed. Also, since the correlations developed 
are empirical, extrapolation to any conditions outside the ranges tested is dan- 
gerous. And finally, the attempts to utilize the correlations in rocket combus- 
tor codes have not been successful. All of the major rocket combustor codes in 
use at Rocketdyne (i.e., TPP, SDER, CICM) have arbitrary multipliers of the 
initial droplet sizes, either as a part of the code or as an input, in order to 
force agreement between the codes and hot-fire engine test data. 

Properties Correlations for Like Doublets 

In addition to all the like doublet "lessons learned" discussed above that are 
applicable in general to rocket engine-type injectors, these studies provide the 
only known corrections or correlations by which we may relate real propellant 
atomization to that of the simulants used in atomization experiments. Although 
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many of the droplet size correlating equations contain liquid properties effects, 
probably the two properties correlations quoted most generally are 
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Ingebo (Ref. 13) 


and 


0 .« y * 333 o' 5 Wolfe and Anderson (Ref. 37) 

jU L l 

Ingebo* s correlation comes from a droplet size correlation equation defining the 
droplet sizes produced by the breakup of a single stream injected transversely 
into an airstream. Wolfe's and Anderson's correlation is based on the breakup of 
already formed droplets in gas streams (i.e., secondary atomization). The 
applicability of either of these relationships to like doublets can be ques 
tioned. In addition, no attempts to establish the effect of liquid properties; on 
droplet size distributions were found in the literature. Also, properties cor- 
relations for unlike doublets, triplets, pentads, or coaxial injectors, or any 
gas/liquid injector apparently do not exist. 

Another aspect of the liquid properties correlations problem that often is over 
looked is the actual values of the properties of the real propellants and the 
simulants at their injection conditions. Since liquids are generally, relatively 
incompressible, since viscosity usually is not considered to be a function of 
pressure, and since density, viscosity and surface tension data for many propel 
lants and test fluids is readily available only at room temperature and one 
atmosphere or at the liquid's normal boiling point (for cryogenics), these room 
temperature and one atmosphere or NBP properties data often are utilized. This 
can cause considerable error. Liquid oxygen is a propellant of considerable 
interest which serves as a good example. For LOX at 134 K (the SSME preburner 
LOX injection temperature) the density increases by 11% and the viscosity 
increases by 52% between 17 and 340 atm (data from NBS Table TN 384) . LOX 
properties are, of course, a fairly strong function of temperature, and choosing 
the wrong temperature (e.g., using NBP data) also can cause great errors. 
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Surface tension is a particularly difficult property for which to find nonroom 
temperature and one atmosphere or non-NBP data. Surface tension is a stron 6 
function of temperature and techniques are available to compute the effect of 
temperature. For LOX, the surface tension changes from 13.2 dynes/cm at its NBP 
of 90 K to 6.4 dynes/cm at its SSME injection temperature of 134 K. As part of 

an attempt to determine the effect of pressure on surface tension, papers were 

found that indicated a very strong effect (e.g., O.K. Rice, "The Effect of Pres 
sure on Surface Tension,: Journal of Chem, _P jiysics, Volume 15, #5, May 1947). 

However, based upon discussions with Prof. A. Adamson and Dr. R. Massoudi of the 
University of Southern California’s Chemistry Department, this effect apparently 
is not due to pressure, but rather to the absorption of gases into the liquid. 

The effect of pressure alone on surface tension should be on the order of a 1% 

increase per 100 atm pressure. This absorption of gases also probably would have 
a great effect on other properties. Since the time available for absorption, 
i e., the time between injection and atomization is so short, very little absorp- 
tion* would be expected. If this is the case, the effect of pressure on surface 
tension should be of little concern. 


Triplet Correlations 

Very little data was found regarding the atomization characteristics for trip- 
lets. This data is synopsized in Table 2. All these investigations were per- 
formed at Rocketdyne utilizing the hot-wax technique. In all these tests, the 
wax was injected into "still" air at ambient pressure. 


As a small part of Zajac’s earlier investigation (Ref. 58), a particular liquid/ 
liquid triplet having all three holes the same size was subjected to atomization 
testing, in order to separately evaluate the droplet size produced by the inner 
and outer streams, wax and hot water were employed. The wax was injected through 
the inner orifice and the water through the outer orifice, and the liquids then 
were reversed on a subsequent test. The only variables investigated were the 
liquids’ velocities, and these were varied in such a way as to maintain a con 
stant mixture ratio. Most of these tests were performed under laminar flow con 
ditions. At high velocity (turbulent flow), the data begins to markedly deviate 
from the correlating equation presented in Table 2. 
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TABLE 2. DROPLET SIZE CORRELATIONS FOR TRIPLETS 
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Figure 6. Correlation of Relative Droplet Diameter with Weber 








Pentad Correlations 


The state of the art for pentad atomization knowledge is essentially the same as 
for triplets. What little data is available was obtained from Rocketdyne wax 
tests. All of these tests were performed by injection of the propellant simul- 
ants into "still" air. The data is synopsized in Table 3. 

As a part of Dickerson’s (Ref. 47) investigation of injector atomization charac- 
teristics, a number of tests were performed on a set of pentad injectors. Drop- 
let size correlating equations were developed relating the mass median drop size 
to the orifice diameters and injection velocities. Separate equations were 
obtained for the inner and outer orifices. Wax and hot water were used as the 
test liquids in a manner similar to the previously discussed triplet tests. In 
addition, droplet size distribution data were obtained. Normalized volume dis 
tribution plots from this work were presented earlier (Fig. 4), and show the 
different distributions obtained for the center and outer orifices. In addition, 
the droplet size distribution equations for this data are presented in Ref. 52. 
As previously discussed, the droplet size correlating and distribution equations 
presented in Ref. 47 are incorrect, and the equations in this latter paper (Ref. 
52) are correct. Dickerson also notes that the quality of the wax spheres was 
poorer than usual for these pentad tests. 

Zajac (Ref. 58) performed a few similar tests and found that the very few higher 
velocity tests were in crude agreement with the correlations of Dickerson. Host 
of Zajac’s tests were at lower velocities and were in great disagreement with 
Dickerson’s correlating equation (Dickerson did not perform tests at these lower 
velocities). The deviation at the low velocities is speculated to be due to vel- 
ocity profile and/or laminar flow effects. Zajac speculates that the flow regime 
of the outer streams is more important than that of the inner stream. 

As a part of the investigation of Mehegan, et. al. (Ref. 55) of gas/liquid injec- 
tors, atomization characteristics were determined for a set of pentads. These 
experiments employed wax and hot gas, with the central orifice always flowing the 
gas. These tests were performed at atmospheric pressure with variations in gas 
and liquid velocity and orifice sizes. No correlating equations were developed. 
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TABLE 3. DROPLET SIZE CORRELATIONS FOR PENTADS 
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TABLE 4. DROPLET SIZE CORRELATIONS FOR COAXIAL INJ 
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Figure 7. Mixture Ratio Influence on Dropsize for Coaxial Injector* 



Investigations by Burick (Ref. 64 and 67) and McHale and Nurick (Ref. 74) util- 
ized hot wax and nitrogen injected into a pressurized chamber to examine coaxial 
element atomization. The combined work of these two studies, along with the pro 
viously discussed work of Mehegan, indicate that recess only reduces D at low 
pressures and/or for large elements. Burick correlated his data as shown in 
Fig. 8. Again, the normalization and logarithmic plotting of the data masks the 
"spread" of the data. Although McHale and Nurick were investigating primarily 
the atomization characteristics of noncircular orifices, they did perform limited 
tests on circular orifices. Their data indicates that increased annulus gas 
dynamic pressure (p V*) reduces droplet size, especially at low liquid 
velocity. References 67 and 74 present droplet size distribution plots. 
Although McHale and Nurick state that recess is a major factor influencing drop 
let size, this conclusion is based upon tests of all of their injectors, which 
are primarily noncircular. The limited testing performed with circular coaxial 
elements indicates a 10 to 20% reduction in drop size as recess is increased to 

R - d . Even the noncircular elements do not show an effect of recess anywhere 
L# 

near as significant as- that found by Mehegan, et. al. 

Falk (Ref. 78) investigated the atomization characteristics of coaxial elements 
injecting wax and hot nitrogen cocurrently (i.e., axially) into a duct flowing 
hot nitrogen. This work utilized the same test apparatus and techniques as the 
analagous work of Zajac (Ref. 70 and 71) on like doublets. One potentially very 

important finding of this work was that the droplet size distribution of these 

coaxial injectors could be described by the distribution function defined by 
Zajac for like doublets. Also, the mass median droplet sizes observed in both of 
these investigations were essentially the same at high relative gas velocity. 
This would seem to indicate that the manner in which the liquid is broken up 
(i.e., the type of injector) has no effect on the ultimate droplet size in the 
presence of a sufficiently accelerating combustion gas. If this is truly the 

case, it is a most important discovery that will direct the course of future 

studies of rocket engine injectors. 
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Some of Falk’s results, showing the effect of simulated combustion velocity on 
droplet size, are presented in Fig. 9. This data indicates that injectors, which 
form larger droplets when no combustion gas motion is simulated (i.e., when 
V =0), show more effect of this gas motion than injectors producing smaller 
droplets. Recognizing this important influence of the relative combustion gas 
velocity on the droplet size, Falk correlated the data in a manner shown in 
Fig. 10. This correlation is based only upon the relative, simulated, combustion 
gas velocity and D q , the mass median droplet size produced by an injector in 
the absence of this gas flow. 

ATOMIZATION SURVEY - FINDINGS, CONCLUSIONS, AND RECOMMENDATIONS 

The state of the art regarding our knowledge of atomization processes is gener- 
ally quite poor. The physics is poorly and, at best, only qualitatively under 
stood. Only very rudimentary quantitative theories exist. The available data 
and correlations are generally of questionable validity and/or utility. Many of 
the most critical parameters are unknown (e.g., combustion gas velocity field, 
multiple element effects) and/or are not simulated in tests (e.g., gas densities, 
real propellant fluid properties, combustion gas motion). This sad state of 
affairs appears to be attributable to two primary causes: the great complexity 

of atomization processes, and the inaccuracies, errors, and limitations associ- 
ated with droplet size measurement techniques. Nevertheless, the available data 
does provide information regarding the importance and relative effects of a num- 
ber of variables on droplet size. 

Probably the most critical of these parameters affecting droplet size is the com 
bustion gas velocity field. This is unfortunate since the actual velocity field 
in a rocket combustor, and in atomization experiments, is unknown. Combustion 
gas velocity also is the one parameter that greatly increases the complexity of 
the atomization assessment problem. This is due to the fact that atomization is 
highly dependent on the combustion gas velocity field, and in turn, the combus 
tion gas velocity field is established by the rate of combustion, which is deter- 
mined by the rate of propellant evaporation, which is highly dependent on how 
well the propellants are atomized (i.e., initial droplet sizes) and mixed. Thus, 
all of these problems are coupled and the solution of any one requires at least 
an approximate solution of each of them. 
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igure 10. Correlation of Droplet Breakup (Atomization) 
Data - Coaxial Injectors (Ref. 78) 



All of the droplet size measurement techniques applied to atomization studies 
have serious limitations and potential and/or known sources of error. Imaging 
techniques measure the spatial concentrations of the various size droplets. Such 
spatial concentrations can be utilized only rarely to define the actual droplet 
size distribution or representative droplet size characterizing all of the drop- 
lets produced by a given spray (temporal distribution). Spatial and temporal 
distributions are often quite different. Thus, the photographic techniques and 
the droplet freezing (i.e., hot wax) technique do not measure the same thing. 

In order to utilize cold- flow atomization data, it is necessary to be able to 
account for the effects of the different liquids’ properties on the droplet 
sizes. The only data available for this purpose applies to like doublets, is of 
questionable validity and applicability, and differs from one investigation to 
another. No methods have been proposed to accomplish this properties effects 
correlation for any gas/liquid injector or for any liquid/liquid injector except 
like doublets. No attempts have been amde to assess injected fluids properties 
effects on droplet size distributions. 

Very little information could be found regarding the atomization characteristics 
of triplet, pentad, and coaxial injectors. Such data, as is available, is pre 
sented along with a representative sampling of the data for like doublets. 

The following actions are recommended for the purpose of (1) improving our knowl- 
edge of atomization processes, (2) developing the droplet size data required by 
the combustor analysis codes, and (3) utilizing the data in such codes. These 
actions are divided into near and long-term approaches. 

Near term: For the most immediate future, it is recommended that droplet size 

data for combustor analysis be determined in the following manner. First, the 
existing data can be utilized (it should be verified first, however) and/or tests 
can be performed to better define D o * the droplet size produced in the absence 
of any simulated combustion gas motion. This can be considered primary atomiza 
tion. Then the data and correlations of Falk (Ref. 78) and Zajac (Ref. 70 and 
71) can be employed to estimate the effect of gas velocity on droplet size. In 
order to do this it is, of course, necessary to estimate the combustion gas 
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velocity where the injector is to be employed. This can be accomplished through 
the use of combustor performance computer codes that generally compute the axial 
velocity of the gas. Thus, at least the major (hopefully) gas phase velocity 
component will be estimated. Since the computed axial gas phase velocity will 
depend on initial droplet sizes, a few iterations of this process may be neces 
sary. That is, the codes can be used to predict V^, which then can be used to 
estimate 5, which will be input to the codes to predict a new V^, etc. 

Another problem in the use of cold flow droplet size data is that it is necessary 
to convert from the test fluids to the real propellants. With great reservation, 
and only because no better information is available, the properties effects cor 
relations of Ingebo (Ref. 13) or Wolfe and Anderson (Ref. 37) are recommended for 
this purpose, when liquid, like impinging elements, are being considered. 


The method described above provides a rudimentary technique for estimating a 
representative droplet size. Drop size distributions in general, and represents 
tive droplet size information for gas/liquid injectors, cannot be estimated via 
this technique due to the lack of data regarding combustion gas velocity effects 
and fluid properties’ effects on atomization. Even when applied to the case of 
like doublets, which have been most extensively studied, this technique may be 
little better than a consistent guessing method. 


In order to better utilize this technique and improve its accuracy the following 
are recommended: 

1. Experiments to investigate gas velocity effects on droplet sizes 

2. Additional tests to better define D q for the injectors of greatest 

interest, especially gas/liquid injectors. Most of the geometric and 

operational variables have not been tested 

3. Experiments to establish fluid properties effects for all types of 
injectors, like and unlike liquid and gas/liquid injectors, and separate 
effects for primary and secondary atomization 
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Such studies and experiments will provide the basis for improvements to atomiza 
tion assessment methods and will establish the nature, feasibility, and desire- 
ability of pursuing the long-term approach. 

Long term: As previously discussed, due to the importance of the combustion gas 

motion on atomization, the problem becomes coupled with those of droplet evapora- 
tion, combustion, and three-dimensional fluid mechanics with momentum and mass 
sources and sinks. Unless some simplifying assumptions are identified earlier, 
the only available solution would consist of a coupling and solution of all the 
equations governing these processes. This would probably involve a long-term 
effort consisting of a number of programs to model (probably with a computer 
code) various parts of the problem, experimentally verify these models, and com- 
bine them in one comprehensive model. Such an approach offers the greatest 
potential for a comprehensive, accurate, proven solution to the problem of spray 
definition for rocket engine injectors. If a satisfactory measurement technique 
exists, experiments with operating, small-scale rocket combustors should be per- 
formed to validate the atomization model. In its ultimate form, such a model 
would include multiple element effects and would predict mixing efficiencies. 
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NOMENCLATURE 


Some atomization nomenclature is defined in text. 

2 

A area (cm ) 

d injector orifice diameter (cm) 

D droplet diameter (microns) 

mass median diameter 

droplet mass median diameter observed when V = 0 

D droplet mass median diameter observed when V ^ - V 

c ° 

f droplet distribution function, (drops/micron) 

f = lim n 
AD-»o AD 

g length over which gas is accelerated in atomization 

accelerating gas flows (cm) 

n number of droplets counted in a giver size range 

n total number of droplets counted 

N cumulative number distribution, 

D 

N (D) = I n 
0 

N flowrate of droplets of size group i 

i 

AP injector orifice pressure drop (Pascals) 

R normalized cumulative volume distribution, v - ^ ^ tot 

W mass rate of flow (kg/s) 

V velocity (m/s) 

V maximum gas velocity 
gm 

V cumulative volume distribution, 

D 

V (D) = E v 
0 

is 

v volume of all drops in a given si ze range ^ in / 

V total volume of all droplets counted 
tot 

V. velocity of drops of size group t 


studies in 
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V simulated combustion gas maximum velocity (gas/liquid injectors only) 
C B 

Y annulus gap for coaxial injectors (cm) 

Z axial spatial coordinate 

a impingement angle 

p viscosity (cP) 

p density 

3 

p concentration of drops of size group i (drops/cm ) 

1 

a surface tension (dynes/cm) 

Subscripts 

g gas (either local chamber gas or injected gas) 

i size group of droplets 

j jet or orifice 

L large drops 

L liquid 

s small drops 
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MIXING 


Cold flow mixing tests frequently have proven to be a significant aid in 

predicting potential performance, or diagnosing problems with rocket engine 

injector components. Cold flow tests are not sufficiently reliable so as to 
serve as a replacement for hot-fire testing, but should be considered as comple- 
mentary to hot-fire tests, aiding in minimizing the number of hot-fire tests 

required to obtain an optimum configuration. In almost every case, an injector 
or element that performs poorly in cold flow testing will not perform well m 

hot-fire testing. However, the counter side of this statement cannot be applied 
universally. An element can be excellent in cold flow mixing, but the combustion 
reaction may override the hydromechnica I mixing provided by the injection 
streams. This effect is most notable with storable hypergolic propellants, when 
a phenomena of reactive demixing "blowapart" is I requently a significant factor 
in combustion performance. There have been other reports of combustion systems 
suffering from reactive demixing, but none have been as well documented as the 
hypergolic reaction systems. 


Aerodynamic forces in the combustion zone also are factors that cannot be simu- 
lated in cold flow mixing tests. Cas forces in recirculation can be strong fac- 
tors influencing mixing and atomization. There are, however, useful correlations 
between cold flow mixing and combustion results, and the relative cost factor 
between cold flow and hot-fire tests generally is a rational reason for utiliz- 
ing cold flow tests as an injector design and development tool. 


The key objective, to establish correlations between cold flow mixing data and 
hot-fire results, requires a large empirical data base as well as a consistent 
assessment of the data and an applied scientific evaluation of the resultant cor- 
relating parameters. Therefore, an assessment criteria was established, which 
allowed compilation of existing cold flow experimental data acquired within the 
industry on element types suitable for LOX/hydrocarbon injector advancement . 


The triplet, pentad, and coaxial element injection devices were selected for 
study based on available hotfire and cold flow experience with LOX/hydrocarbon 
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propellants. The cold flow correlating parameters used for these devices were 
identified and an extensive literature survey conducted to obtain related cold 
flow data. Data from the literature search was compiled into a displayable for- 
mat. The information then was plotted by the appropriate correlating para- 
meter(s) against mixing efficiency, a standard measure of cold flow performance. 

In addition to the literature survey, five impinging triplet elements, one pentad 
element, and three coaxial elements were fabricated for cold flow testing. The 
sizing of these elements encompassed designs for both preburner (gas generator) 
and main injector mixture ratios at high chamber pressure. The propellant com- 
binations were LOX/methane (gas/liquid), LOX/RP— 1 (liquid/liquid), and LOX/pro- 
pane (liquid/liquid and gas/liquid). The low-pressure cold flow mixing test pro- 
gram was conducted with these elements at several flow conditions. Measures of 
mixing efficiency were established and plotted as a function of mixing param- 
eters. Maps depicting mixture ratio-normalized mass distribution were con- 
structed from the cold flow tests to provide a good visual indication of relative 
mass and mixture ratio concentrations for the different element types. 

INJECTOR MIXING CORRELATING PARAMETERS 

Mixing correlation parameters are mathematical expressions based upon injector 
element geometry and flow conditions. Their utility as injector design criteria 
depends upon (1) their ability to be related to mixing efficiency and (2) the 
existence of optimum values of these correlation parameters at which mixing will 
be maximized. 


Numerous correlating parameters have been proposed for different injector config- 
urations, propellant conditions, and hot fire related operating conditions. The 
scientific basis for the parameters generally has been derived from momentum and 
stream diameter relationships of the injection element. A survey of available 
literature showed that of these relationships, most correlating parameters were 
derived for liquid/liquid impinging- type injectors. Many of the experimenters 
have established formulas to plot data from numerous test conditions on a single 
curve, or at least, within a family of curves. 
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The correlating parameters used in the literature survey data reduction and in 
the subsequent low-pressure mixing tests are presented in Table 5. A descrip- 
tion of these important parameters is discussed below. Illustrations of the 
three element types studied under this program (coaxial concentric tube, triplet, 
and pentad) are presented in Fig. 11 through 12, with the appropriate terminology 
and physical parameters identified. 

Rupe Factor/Rupe Number 

The best example of an injector correlating parameter for mixing criteria is the 
Rupe Factor, or Rupe Number, developed for use on unlike impinging doublets ele- 
ments. This basic expression (F.q. 2) primarily was developed in the '50s by its 
namesake, Jack Rupe of Jet Propulsion Lab (JPL). He ran a great number of cold 
flow mixing tests and conducted related hot fire experiments. Applying stream 
momentum and diameter ratios, he developed an expression, since referred to as 
the "Rupe Factor," which indicated the best mixing when it equalled unity. This 
parameter also can be expressed as the diameter ratio over the momentum ratio. 
Since this expression is a ratio, the mathematical range of this factor from ,.cto 
to one is the same as from one to infinity, which is difficult to interpret. lot 
this reason, the expression has been revised to the "Rape Number" (Eq. 3), which 
has a total range from zero to one and an optimum -a Luc of 0.5. 

This expression has been utilized widely lor sizing of unlike doublets and has 
demonstrated good correlation over a vide range of conditions. This does not 
mean Chat a Rupe number of 0.53 reflects a certain quantitative level of mixing 
efficiency, but that in sizing an element for a given design, mixing, for most 
cases, optimizes very near the 0.5 value. 


Momentum Ratio 


Other element types have been analyzed in a similar manner as the Rupe Number, 
and modified momentum/diameter relationship expressions have been derived for 
triplet and pentad impinging element patterns. These parameters are based on 
more limited cold flow data and virtually no hot. fire data, and should he used 
more cautiously in universal application than the doublet expressions. 
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OUTER TUBE I.D. (FUEL SLEEVE) 



Figure 12. Impinging Triplet and Pentad Elements 


As previously mentioned, the primary root for almost all impinging element mixing 
parameters is the momentum ratio. As a general rule, the momentum ratio always 
is expressed as the oxidizer total momentum over the fuel total momentum regard- 
less of the number or placement of oxidizer streams relative to fuel streams 
within the element. Relating this ratio to the values available to the designer, 
we have the form of momentum ratio as shown in Table 5, (Eq. 1). There is no 
design optimum for this parameter and, again, this is a ratio with tneoret leal 
values from zero to infinity, where values over one indicate that the oxidizer 
lias higher momentum than the fuel. 


Klverum-Morey Factor 

The equivalent of the Rupe Factor for triplet and pentad elements was developed 
by Rupe’s colleages, Elverum and Morey, and is Hised also on momentum/diameter 
(area) relationships as shown in Table •>, (Eq. 4). For the triplet element, with 
two outer angled streams and a central axial stream, the relationships are set Uo 
inner and outer streams rather than oxidizer and fuel streams, since both 
fuel -oxidizer-fuel and oxidizer-fuel -oxid izer triplets are in general use. 


For liquid/liquid triplets, within the range of study by Elverum and Morey, the 
optimum value for this expression was 0.6b. Triplet injectors have been used 
most commonly Cor hypergolir. storable propellants, and use of the ELverum-Morey 
Factor has been successful under these conditions. For the nominal mixture ratio 
of Liquid oxygen and liquid hydrocarbons*, very little data has been available. 
A modified Elverum-Morey expression. Table 5 ( Eq . 5), was designed for pentads 

and has a purported optimum value of 2.75. 


*The typical mixture ratio for storable propellant, combinations, such as NTO/MMH 
or UDMH/IRFNA, is between 1.5 and 2.5 ox/fu for main injector operation. The 
mixture ratio for liquid/liquid LOX/hydroc arbon propellants, i.e., RP-1/L0X, is 

optimum near 2.8 ox/fu for main injector maximum Isp, and near 0.4 ox/fu for 
fuel-rich preburner (turbine drive combustor) applications. 
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Pe ne t ration Factor 


This parameter has been developed for gas/liquid triplet injectors where two 
liquid streams impinge on a central gas stream. It relates the predicted 
penetration of the liquid streams to the central gas flow. Optimum mixing is 
predicted if the liquids barely penetrate to the center, with liquid droplets 
being sheared off and entrained by the gas flow on the way. The penetration fac- 
tor is presented in Table 5, (Eq. 6). 

A value of 0.5 is the theoretical optimum. Lower numbers infer that the liquid 
is being deflected away by the gas or is not fully penetrating the gas stream. 
Over penetration, on the other hand, produces a liquid fan within the gas flow, 
which also reduces the uniformity of gas/liquid mixing. This factor was created 
from a combination of analysis and cold flow experiments, and hot fire data 
appears to support the basic premise. Pentads and other impinging patterns with 
liquid streams impinging on a central gas core also could be expected to corre- 
late with some form of the penetration parameter. However, data is limited for 
these applications . 

The use of this factor for the reverse case of gas streams impinging on a central 
liquid, or any other extremes in the density relationships, is questionable. 
Triplets with the gaseous reactant in the outer streams have been used in num- 
erous cases, but there is little data on any correlating parameters. Some 
limited information suggests that high levels of gas to liquid momentum ratio are 
beneficial to the mixing process in impinging element injectors. 

Velocity Head Ratio 

Another parameter that does not have a stated optimum value is the velocity head 
ratio shown in Table 5 (Eq. 7). This roughly relates to the very practical 
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consideration of "Delta P" ratio, or pressure drop ratio. The usual starting 
point in an injector design is based on the desired level of pressure drop at the 
design flowrates. Isolation between feed system and chamber pressure distur- 
bances generally dictates a desire for a high level of injection orifice pressure 
drop, and system pressure limitations would like a low pressure drop. A compro- 
mise solution usually results in an injector delta P of about 15 to 20 % of cham- 
ber pressure, and an initial starting point would he for both oxidizer and fuel 
systems to be roughly the same value. Therefore, an injector design that has 
velocity head ratios significantly distant from i would require some compensa- 
tion in design approach (i.e., supplementary orifices, etc.). 

As mentioned previously, there is no theoretical optimum for the velocity head 
ratio, but the values close Lo 1 are desirable for system integration. Many 
times, sizing the injection orifices to optimize on; of the other parameters will 
result in an unacceptable level of velocity head ratio. For this reason, the 
velocity head ratio should be computed at the sam< time as the other parameters, 
and evaluated and adjusted concurrently. 

Coaxial Parameter 


The gas/liquid coaxial concentric tube injector element has had wide, success! u 1 
usage for hydrogen/oxygen combustion. Cold flow and hot fire experience with 
this element still has not provided a good correlating parameter. In this ele- 
ment, typical design practice has been to provide a low-velocity central liquid 
Stream (liquid oxygen) sheathed by a high-velocit > gas flow (gaseous hydrogen or 
fuel-rich preburner gases) as shown in Pig. 11. Mixing and atomization are pro- 
vided primarily by the shear forces between gas and liquid and by the momentum of 
the expanding gases. 


Recessing the liquid stream upstream of the exit "lane of the outer (gas) stieam 
is popularly held to increase both atomization and mixing. Cold flow testing has 
not established a strong correlation with this practice, although hot lire 
results generally reflect a performance increase that usually is accompanied oy 
an increase in face heating.. 
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Increasing the gas velocity (relative to the liquid velocity) generally improves 
mixing. This design approach should not be employed blindly, since some refer- 
ences suggest that mixing can be impacted adversely by velocity ratios that are 
too high. This would tend to suggest that some correlating parameter for optimi- 
zation may be possible. Very high gas velocity apparently can reverse the gas 
liquid relationship, "blowing out" the center of the spray and dispersing excess 
liquid to the outside of the spray cone. 

A review of existing data, as a part of this effort, indicates trends that may be 
useful for providing a general optimizing expression *' or the coax element. Falk 
and Nurick of Rocketdyne (NASA CR-72703 R-8361) have suggested the coaxial para- 
meter presented in Table 5, (Eq. 8). However, no optimum value of this parameter 
has been established. One of the objectives of the remainder of this program is 
the establishment of a coaxial element mixing parameter. 


MIXING TEST METHODS 

Liquid /Liquid Mixing Test Methods 

The liquid/liquid testing for mixing efficiency is relatively easy and low cost, 
if facilities are available. The procedure for liquid/liquid mixing utilizes a 
grid-like sample device, which ducts the individual position captured liquid into 
an appropriate sample container (Fig. 13a). This technique utilizes two 
immiscible liquids as propellant simulants, typically water and a high-density, 
low-vapor pressure solvent such as 1,1,1-trichloroethane. The fluids collected 
in the sample tubes separate by the variation of density and their quantities in 
each tube are measured (Fig. 13b). Typically, the sample grid represents hundreds 
of data points, and a computer data reduction process is required to provide 
meaningful quantitative data. 

Different fluid combinations have been employed for liquid/liquid mixing in an 
effort to better match injected reactant conditions, while addressing concerns 
for toxicity, flammability, and general questions oi safety, convenience, and 
cost. Other solvents used for these purpose: have included many of the lower 
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vapor pressure "freon" compounds, perchlorethylene (a dry-cleaning solvent), as 
well as fuel-type hydrocarbon liquids. At least one past program at Rocketdyne 
utilized a water/brine system, with the mixture ratio of the sample determined by 
an electric salinity meter. Data acquisition using this method was significantly 
slower than the immiscible fluid method, and accuracy was poor in the low mass 
flow outer zones. 

Gas-Liquid Mixing Test Methods 

Gas-Liquid mixing tests are significantly more time consuming than the liquid- 
liquid mixing, which probably is the reason that gas-liquid data is more 
limited. A gas-liquid mixing measurement system has been utilized extensively at 
Rocketdyne for hydrogen/liquid oxygen concentric elements (with the gas annulus 
surrounding the liquid core). The schematic of the process is shown in Fig. 14. 
The sample element is installed at the "head end" of a transparent, pressurized 
chamber, with a traversable probe mounted at the desired sampling plane. Water 
typically is used for the liquid oxygen simulant and a nonreactive gas simulates 
the hydrogen fuel (or hydrogen— rich hot gas in a staged combustion cycle). 
Typically, the gas used is nitrogen, sometimes dilute! with helium to provide a 
desired density. Gas density is controlled by tank back pressure, and the mix- 
ture of gases supplied. A "base bleed" gas usually is supplied through the face 
around the injection element to minimize recirculation from the injected flow, 
and to simulate partially the axial gas flow present n a combustion chamber. A 
tracer gas (frequently oxygen) is included in this base bleed flow to allow this 
local gas flow to be measured and extracted mathematically from the measured 
element gas flow in each sample. 

The sample is extracted from the gas-liquid element flowfield by the use of a 
sharp edge probe that can be positioned in the desired sample area. The liquid 
spray in the sample zone is collected physically by the opening in this probe, 
and accumulated in a sample container over a measured time period. The gas flow 
flux in the sample zone is determined from the relationship between total and 
static pressure (corrected for the liquid in the two- phase flow). The gas mea- 
surement may require a second correction for the entrained base bleed flow, and 
the data for this correction is obtained from an "on-line" gas analysis technique. 
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As might be deducted from the preceding description, each sample requires a 
sufficient time to stabilize the required readings and collect the liquid. When 
compared to the hundreds of sample points simultaneously obtained in the 
liquid-liquid testing, the increase in test time for gas-liquid testing is 
readily apparent. Testing with concentric elements permits a reasonable 

assumption of circular symmetry, allowing a reduced number of required sample 
measurements. However, the more complex "fan” shapes of gas-liquid triplets and 
pentads require careful study of the sample locations, and require more sample 
points than for a co-ax test. Previous work with triplets and pentads in a 

gas/fluidized solid system, Ref. 30, and triplets in a gas-gas system have 
indicated the shape of expected mass distribution, and show that numerous sample 
points are required to characterize these element types. 

COLD FLOW MIXING DATA REDUCTION 

The data reduction procedures for the liquid-liquid and the gas-liquid cold flow 
mixing tests are very similar. As in the testing itself, the data reduction for 
the liquid-liquid testing is a bit more straigat forward. The total sample grid 
usually encompasses all the injected flow, and the grid openings usually have no 
open spaces between them. Therefore, the collected totals should equal the 
injected totals, thus providing a good cross-check on the data. This is the 
first factor computed in gas/liquid mixing tests, and is referred to as the 
"collection efficiency." 

Collection Efficiency 

To calculate the collection efficiency of the test system, fluid input values are 
compared with fluid collected values. The input values of mass flow rate W is 
frequently calculated theoretically by the Injector Vressure Drop Equation (9), 
based on previous cold flow resistance calibration of the test model. 

Winput = l Nd2 C d (2 pAP) l/2 (9) 
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where 


W = mass flowrate 

N = total number of oxidizer or fuel holes 

D = diameter of orifice 

Cd - dimensionless discharge coefficient as determined 

from the calibration flow test 


p = density of simulant 
AP = injector pressure drop 


If direct, flow measurement capability exists in the cold flow mixing facility, 
the values from these measurements are used. 


The collected values of mass flowrates are calculated from the test data; sum 
ming all of the individual sample measurements: 

w collected = ^ ( L<n 


where : 



P 

- density 


Q 

- local corrected 

sample voluim 

t 

- collection t ime 

in seconds 


Collection efficiency of the system is calculated then by. 


ft collected ( \ \ ) 

col W input 

where a value of ”1" represents perfect collection efficiency. Large deviations 
in the collection efficiency would indicate problems in the system or the data 
for the testing. Unfortunately, collection efficiency rarely is included in 
reports of mixing tests and in some cases may not even be calculated. Liquid/ 
liquid mixing testing is relatively simple and collection efficiency generally 
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is not needed or obtained. However, the much greater complexity of gas/liquid 
testing requires the "check” on test methods and procedures that collection 
efficiency provides. 


Mixing Efficiency 

The most meaningful expression for assessing mixing efficiency is the E m <E-sub 
m) value proposed several years ago by Jack Rupe at JPL. This is an expression 
for the mass mixture ratio distribution of the samples based purely on the rela- 
tionship of the samples to the overall mixture ratio with no regard to such fac- 
tors as theoretical stoichiometry, etc. 

This value is computed as a mass weighted summation of the mixing errors in all 
the samples. In practice, it is computed as a summation of decrements based on 
how far the mixture ratio of each sample deviates from the overall mixture ratio, 
and weighted by the mass fraction of each of these r.amples. The range of this 
expression is from zero to 100%, with 100% indicating all samples are the same 
mixture ratio, and zero indicating the samples are all one component or the other. 


The nominal form for computation of is expressed by. 


'm 


where 


R 

R sb 

M R sb 

**F S a 

R sa 


100 


1 IMF 


R-R 
sb R 


sb 


R-R 


4 IMF 


sa 


sa R-l 


Efl, = mixing efficiency from 0 to 100% 


overall mixture ratio as expressed by weight flow 
oxidizer/weight flow total 

mixture ratio of sample below overall mixture ratio 
mass fraction of sample below overall mixture ratio 
mass fraction of sample above overal mixture ratio 
mixture ratio of sample above overall mixture ratio 


(12) 
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Each Local sample that Ls not at the overall mix.uo ratio thus provides a mixing 
efficiency decrement proportional to how far it is from the nominal mixture 
ratio, and what mass fraction of the total I low it represents. For example, if a 
sample representing 50% of the total mass has a mixture ratio fraction of 0. ^ 
when the overall is 0.70, the total mixing loss from this sample is 

100 (0.5 - 7 °q~ 7 § ~~ > = loss Ln mixi "8 efficiency 

This factor is much more sensitive to mixing deficiencies than combustion 
efficiency-related factors, which are "rounded off" by theoretical curves and the 
relationship between test mixture ratio and s toich Lornet ric mixture ra^io. 


Mixing Limited C~Star 


A frequently used parameter to describe mixing test results is mixing limited 

C-star or * mix (ETA C-star mix). This can be applied only to tests for a 
c 

specific reactant combination, and actually only for an assumed chamber pres- 
sure. It is a prediction of the expected hot fire C-star efficiency (assuming 
total vaporization). The product of vaporization efficiency and mixing limited 
C-star efficiency is the predicted combustion efficiency. 


At Rocketdyne , the mixing limited C-star is computed by a single stream tube 
performance model technique. The computer program is provided with a theoretic il 
C-star function and the theoretical C-star value -M/sec) is calculated for each 
sample mixture ratio. Each sample collected mass is multiplied by the sample 
C— star, and these products are summed for the entire sample. This answet is 
divided by the total mass collected to provide the mixing limited C-star. 

C* x Mass + C* x Mass 0 + C* x Mass 

l .1 2 2 N M T) 

Mixing Limited C star — Total Sample Mass 
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The mixing limited C-star efficiency then is determined by comparing this value 
to the theoretical C-star for the overall mixture ratio: 


mixing limited C* 
theoretical C* 


(14) 


A C-star efficiency of one indicates that at uniform mixture ratio, ^ = 1, 

mixing limited C-star is equal to theoretical C-star. This parameter is used 
make a rough estimate of performance potential for given operating conditions of 
certain mixture ratio and mixing efficiency. 


REVIEW OF EXISTING DATA 

An extensive literature search was conducted on pas. experience in determining 
and evaluating mixing efficiency for triplet, pent id , and coaxial elements. 
Numerous document references were accessed and reviewed, and a bibliography of 
the pertinent reports reviewed is presented herein. The intent of this search 
primarily was to find reports containing quantitative cold flow mixing test data 

for these injectors. 


The literature search yielded fewer reports than had been anticipated, although 
several valuable references were encountered. The ibundance of data involved 
liquid/liquid impinging doublets followed by liquid/liquid triplets and pentads. 
Gas/liquid reports were almost entirely limited to coaxial elements and presented 
little data regarding gas/liquid triplets or pentads. 


The data from each report was re-reduced in order to provide a uniform basis for 
comparison. In each instance, the objective was to obtain as close to raw data 
as possible from the information in the report. Using a computer program 
designed for this task, a table of injection parameter relating to measured per- 
formance was constructed. Information from each repo *t thus was computed in the 
same consistent fashion for best comparison of results. 


As expected, many important test conditions typically were omitted from the 
reports, such as the distance from the injector fac • to the sample plane, the 
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relative size of the sample grid, and the number of sample points in the test 
plane. For gas/liquid coaxial element data there had been controversy on the use 
of averaged data for sample grid points, and the reports typically did not ela- 
borate on data reduction methods. With these limitations and constraints in 
mind, the data was analyzed and reviewed for some generalized conclusions. 

The data was extracted from all the reports that had usable mixing data and has 
been prepared in summary chart form (Table 6). The data has been organized by 
element types and propellant condition (i.e., gas-liquid triplet, liquid-liquid 
pentad, etc.). All of the normally used injector sizing and operating parameters 
are displayed (if they were available or calculable from the report informa- 
tion). Where a report provided information on more than one element type or pro- 
pellant condition combination, it has been listed in appropriate multiple loca- 
tions in the charts, with cross-reference to the other elements. These charts 
are intended as a summary reference source, rather than a quick graphic compari- 
son, and a review of data comparing similar configurations can be accomplished 
with minimum confusion. Most of the data also is presented elsewhere in this 
report in graphic form, with mixing efficiency plotted against the common 
injection parameters. 

Triplet - Liquid/Liquid 

Two documents for liquid/liquid triplets, were found each containing significant 
single element data on several configurations (Ref. 5 and 7). The data was rela- 
tively consistent and indicated a reasonable correlation with the Elverum Morey 
Factor (Fig. 15 and 16). These plots depict the Elverum Morey Factor on a 
logarithmic scale since this factor is computed as a ratio. In both references, 
it can be stated generally that maximum mixing efficiency occurs near the 0.66 
value for the factor. Elements with near the same orifice diameters appear to 
provide the highest maximum mixing efficiency, and multiple elements reflect 
interelement mixing with higher average values and reduced sensitivity to t he 
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ELVERUM-MORET FACTOR 



characteristics would appear to be warranted since liquid/liquid hydrocarbon 
mixtures favor a reverse triplet configuration at main chamber mixture ratios. 

Triplet - Gas/Liquid 

One report was found on LOX/hydrogen work (Ref. 34) which provided data on a 

liquid/gas/liquid configured element. Suprisin^ly, the penetration factor, 

designed for liquid/gas/liquid elements, did not produce the desired correlation 

of maximum mixing efficiency (E ) at the 0 !> theoretical optimum value 

m 

(Fig. 17). Visual aids from the report depict i he gas/liquid normalized mass 
flux profiles for each of three cold flow tests. Figure 18 depicts represent a 
tive samples of those three tests. The sample mixture ratio is equivalent to the 
overall inlet mixture ratio where the dashed lines (gas) intersect the solid 
lines (liquid). It can be inferred from the distribution plots that the balance 
of gas and liquid was optimum at the under- penetrated condition (penetration 
factor 0.4), which contributed to the maximum- measured mixing efficiency. At 
penetration factors greater than 0.4, the gas blowout produced by the impinging 
liquid jets was visible. This contributed to the poorer mixing efficiency noted. 

The Elverum-Morey criteria for this element, shown in Fig. 19, did reflect a 
correlation between the 0.66 optimum value and the peak mixing efficiency. in 
this test, the oxidizer- to- fuel density ratio was over 600, markedly removed from 
the design application range of 1.7. These parameters bear additional testing 
since there are good designs for liquid/gas/liquid elements in LOX/ hydrocarbon 
gas generators and preburners. 

Pentad - Liquid/Liquid 

Documents obtained with mixing data for liquid/liquid pentad elements consisted 
primarily of reverse configuration* element studies (Ref. 3,. 4, and 10) . In these 

* A reverse pentad generally is considered to have the denser liquid (oxidizer) 
in the outboard streams and the less dense liquid (fuel) as the centrally located 
stream. 
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Depicting El verum- Morey Factor Influence 



studies, the overall level of mixing efficiency was generally good. Single 
element characteristics did not adhere to the Elverum-Morey theoretical optimum 
very well for the large element tests shown in Fig. 20 (Ref. 3), although the 
multi-element tests did show peak mixing efficiency near the 2.75 optimum value 
for the same experimentors . This is either a result of secondary mixing enhance 
ment from the multiple element configuration or is indicative of absolute size 
limitations in parameter application. Other data presented in Fig. 21 and 22 
indicate some small degree of correlation with the 2.75 optimum parameter value. 

Pentad - Gas/Liquid 

The volumetric unbalance realized with gas/liquid propellant combinations fee 
quently dictates the use of pentad (four on one) elements. With the gaseous 
reactant on the four outside elements, this bears some resemblance to an imping 
ing concentric element. 

With the gaseous component of the reaction system in the center stream, the case 
resembles an extension of the liquid gas liquid triplet where a form of the pene 
tration factor becomes the most likely mixing parameter. 

One document was located with gas/liquid pentad data (Ref. 31), which includes 
test data for both configurations, liquid gas liquid and gas-liquid- gar. . This 
data was replotted against three different parameters, momentum ratio, Elverum 
Morey ratio, and penetration factor. 

Both pentad configurations showed improved mixing characteristics with increased 
oxidizer (liquid) momentum (Fig. 23), regardless of the orientation of the oxi 
dizer stream(s). This is not understood fully since prior experience on other 
programs, such as the gas/fluidized- sol id program (Ref. 30), indicated contrary 
relationships, i.e., an increase in performance with a reduction in momentum of 
the central-fluidized stream with maximum performance occurring at a relatively 
high gas to liquid momentum ratio. 
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Extrapolating the liquid-gas-liquid test data along th- Elverum-Morey curve. Fig. 
24 suggests a trend toward the 2.75 optimum value for pentads, whereas the gas 
liquid-gas data do not obey the parameter functions. Extrapolating the penetra 
tion factor data for the gas-liquid-gas element may indicate a trend toward the 
0.5 optimum value (Fig. 25). The liquid-gas-liquid element apparently does not 
adhere to the penetration factor function. 

Concentric Coaxial Element 

Several report references were obtained in the literature search containing cold 
flow mixing data for coaxial elements. Some of these were from the Space Shuttle 
Main Engine (SSME) Program. The mixing data from these sources were plotted 
against the conventional parameters applicable to coaxial injectors, namely LOX 
post recess and velocity ratio. 

In most concentric element configuration, relatively large improvements in mixing 
are anticipated as the central tube (oxidizer post) recess is increased to one 
liquid stream diameter. Data presented in Fig. 26 Ref. 7) depict less effect 
than had been expected. The curve indicates poor overall mixing efficiency (E^ 
= 50 to 65%) with very little improvement obtained as recess is increased. How- 
ever, Falk and Burick report in their studies (Ref 19) that cup recess does 
improve mixing. This conflict needs to be resolved by additional testing, espe 
cially in the areas of hydrocarbon fuels. 

The influence of gas-to-liquid velocity ratio on the level of mixing efficiency 
is presented in Fig. 27 and 28, depicting the characteristics of SSME LOX/ 
hydrogen preburner and main injector elements in cold flow test. In these fig- 
ures, mixing efficiency is consistently high. Prop, llant density matching was 
achieved for these tests, which also resulted in nominal matching of hot fire 
(design range) velocity and momentum ratios simultaneously. 

Additional tests conducted by Rocketdyne (Ref. 31 and 7) are presented in Fig. 29 
and 30, respectively, depicting the effects of velocity ratio on mixing 
efficiency. The latter figure shows the influence of gas-to-liquid density ratio 
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Figure 26. Mixing Efficiency for Coaxial Element (NASA CR-R-9270) 
Depicting LOX Post Recess Influence 
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as well, and clearly indicates that higher gas-to-fuel density ratios produce 
higher mixing efficiencies for a given velocity ratio. This relationship 
strongly suggests that a velocity— density product, such as momentum ratio, will 
not peak at an optimum value, but will approach ideal mixing as the gas momentum 
continuously increases. For this reason, an alternate parameter (Table 5. Eq. 8) 
has been considered in an effort to characterize the data with a single 
expression. The coax parameter (Ref. 19) was applied to the SSME preburner and 
main injector data as shown in Fig. 31 and 32). Because of the high overall 
mixing efficiency of that data, no predominant trends were evident. 


LITERATURE SURVEY CONCLUSIONS AND RECOMMENDATIONS 

As a result of the literature review and data examination, most of the initial 
impressions regarding the state of the cold flow data have been confirmed. Large 
discrepancies exist in test results noted between the various experimenters, and 
there does not appear to be any proven correlating parameters for coaxial element 
mixing efficiency. In general, the available data is insufficient to confidently 
confirm or establish the optimum value of the correlating parameters for imping- 
ing elements. 

Although the gas/liquid triplet element has significant potential for future 
liquid-oxygen/gaseous-hydrocarbon propulsion systems, very little quantitative 
data exists to either support design calculations or provide correlating 
expressions for combustion modeling. Most hydrocarbons considered for advanced 
booster applications will be delivered to the injectors as warm or hot gas with 
densities relatively high as compared to hydrogen or combustion gases used in 
current concentric element injectors. This higher density favors impinging 
elements rather than the concentric element. The gas annulus gap required for 
the denser fuels in a coaxial element injector may approach small absolute values 
that ultimately result in poor concentricity and element contamination problems. 
Greater emphasis should be placed on obtaining mixing data on gas/liquid 
impinging (especially triplet) elements. 
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Figure 31. Mixing Efficiency for Coaxial Elements ( 1L PT 73 37) 
Depicting Falk and Burick Parameter Influence 
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COLD FLOW MIXING TESTS 


A series of single element, cold flow, mixing tests were performed to (1) estab- 
lish the validity of the correlating parameters as predictors of mixing effi- 
ciency, and (2) to define the mixing characteristics of particular LOX/hydrocarbon 
injector designs. These tests were performed in the Rocketdyne Engineering Devel- 
opment Laboratory and are described in detail herein. 

Injector Elements 

Test hardware for liquid/liquid and gas/liquid cold flow testing was designed and 
fabricated at Rocketdyne under Task II program objectives. A set of manifolds and 
replaceable single-element injector inserts were fabricated per Fig. A-l 
(Appendix A). The physical data for each element, as well as the applications, 
operating ranges, and design rationale are presented in Table 7. 

Specific selections and sizing of preburner element types were based on a subscale 
hot-fire evaluation program conducted by Rocketdyne under NASA contract NAS8-33243 
(Ref. 32). A wide range of 2— inch diameter LOX/hydrocarbon preburner injectors 
were fabricated and hot-fire tested, which included (1) a fuel-rich LOX/methane 
coaxial element, (2) a fuel-rich LOX/RP-1 triplet, (3) a fuel-rich LOX/methane 
triplet, (4) a fuel-rich LOX/methane pentad, and (5) an oxidizer-rich LOX/methane 
pentad. Based upon those tests, several larger subscale preburners were selected 
for fabrication and delivered to NASA. 

Main injector gas/liquid coaxial elements for LOX/methane and LOX/propane pro- 
pellant combinations were derived from design criteria or analyses performed by 
Rocketdyne under NASA contract NAS8-33206. Under this program (Ref. 33), a high- 
performance LOX/methane main injector was designed, fabricated, and delivered to 
NASA for hot-fire evaluation. 

The electro-deposited Nickel (EDNi) preburner triplet pattern was chosen because 
of Rocketdyne's on-going independent hot-fire research utilizing compact, high- 
element density injectors. The rationale for micro-orifice patterns is to achieve 
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high performance In short distances and realize improved combustion gas spatial 
temperature uniformity across the preburner discharge. 

The two liquid/liquid injector triplet elements were configured to obtain an 
Elverum-Morey factor of 0.66. The closeness of orifice diameter sizes at an 
Elverum-Morey of 0.66 products a generally high level of mixing efficiency. 

Mixing Test Equipment. Procedures, and Data Reduction 

The cold flow mixing tests were performed in the Atomization and Mixing tests 
facility of the Engineering Development Laboratory at Rocketdyne. The test appa- 
ratus consists of two separate units, one for testing liquid/liquid elements and 
one for testing gas/liquid elements. These apparatus and the associated test pro- 
cedures are each described below for each type of test. The more general discus- 
sion of apparatus and procedures contained in a preceeding section of this report 
also applies to these tests. The data reduction and compilation techniques and 
software are presented in Appendix B. 

Liquid-Liquid Mixing . The liquid-liquid mixing test facility utilizes a sample 
collection system consisting of a 13 by 20 sample grid of 0.318 cm square tubes as 
shown in Fig. 13a. This grid size thus provides resolution of the spray into 260 
regions, and is a good compromise between the desired small sample size, and the 
practical flowrates in small tubes. This small grid size is designed espe- cially 
for single element tests. The sampling time is controlled by air actuated 
shutters that divert the flow away from this grid before and after the desired 
sampling period. Each of the samples is ducted to 50 milliliter graduated 
cylinders where the two liquids separate. The oxidizer simulant is colored to 
ensure readability (Fig. 13b). The quantities of each fluid are measured and sub- 
sequently input to the computer for data reduction. 

Propellant simulant liquids used for these tests were 1-1-1 trichloroethane for 
the oxidizer simulant and water as the fuel simulant, providing a density ratio of: 

Density 1, 1. 1-trichloroethane 82.6 

Density water * 62.4 = 1,32 

which is not too far from L0X/RP1 value of 1.42. 
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LIQUID-LIQUID MIXING TESTS 


Test Results 


Three different simulated injection elements were tested in the liquid-liquid 
mixing program. All of these elements were triplet configuration single-element 
models. Two elements (No. 1 and 6) were designed for preburner, or gas generator 
flowrates with "fuel-rich” design mixture ratios for low turbine inlet tempera 
tures. Injector No. 1 was designed for a liquid oxygen/liquid RP1 gas generator, 
and injector No. 6 was designed for liquid oxygen/liquid propane reactants. The 
other element (injector No. 8) was designed as a main chamber element for a liquid 
oxygen/liquid RP1 reactant system. Detail dimensions of these elements are shown 
on Table 7. 

A summary of the liquid/liquid mixing test conditions, correlating parameter val- 
ues, and measured mixing efficiency is presented in Table 8. The nominal mixture 
ratio (NOM MR in the table) is the simulant's mixture ratio that was obtained when 
the simulant's injector pressure was set to provide a flow that would approxi- 
mately match the momentum ratio of the real fluids at their design operating con 
dition. Variations of approximately 20% on this mixture ratio also were tested. 
An attempt also was made to determine the effect of the grid location (i.e., col 
lection distance) on measured mixing. The results for each element are discussed 
in detail in the following. 

Triplet No. 1 . Triplet No. 1 is a single element model of a fuel-rich gas gener 
ator injector for liquid oxygen and liquid RP-1, a kerosene-based liquid hydro 
carbon fuel. This is a conventional triplet configuration with two fuel streams 
impinging on a central oxidizer stream. Fuel orifice diameters are .14 cm and 
they are angled inward at the traditional 30 degrees from axial. The oxidizer 
diameter is .116 cm and the centerline impingement distance is .64 cm from the 
injector face (further specifications are shown in Table 7). 

Triplet No. 1 exhibited disappointing performance in the cold flow mixing tests. 
Mapping the mass distribution from the test rerults indicated significant 
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maldistribution and nonsymmetry, indicating either a manufacturing or manifolding 
defect. The data should not be considered representative of triplets of this con- 
figuration, at these operating conditions. The data does demonstrate how cold 
flow measurement techniques can identify and quantify fabrication errors that may 
not be detected through normal quality control techniques. 


One of the most arbitrary aspects of mixing assessment tests is the choice of the 
collection distance. In order to assess the importance of this choice, triplet 
No. 1 also was tested at a collection distance of 2 5 cm, that is, at half the 
previous collection distance. Flowrates were approximately the same in both 
tests. A 20% reduction in mixing efficiency was observed at the smaller collec- 
tion distance. This one test of collection distance is insufficient to form any 
conclusions. However, it does indicate a need for further investigation. 

Triplet No. 6 . Triplet No. 6 is an element from a &; s generator/preburner design 
for fuel rich operation with liquid oxygen and liquid propane. It is a single 
element model of a conventional triplet configuration utilizing two outer fuel 
streams impinging on a central oxidizer stream. Turn fuel orifice diameter is 
.2 cm, and the impingement angle is the traditional 30 degrees each side. The 
oxidizer orifice is .13 cm diameter and the centerline impingement distance is 

.64 cm. 

The cold flow mixing test of this element indicated b tter symmetry than for trip 
let No. 1, although the mixing efficiency value was :till relatively low (in the 
60 to 70% range). The low mixing efficiency undoubtedly reflected the relatively 
large mismatch between oxidizer and fuel orifice diameters. The oxidizer orifice 
is too small relative to the two fuel holes (each on. over double the area of the 
oxidizer orifice). Some minor skewing was evident in the mass distribution pat- 
tern, but it was felt that it was not particularly significant. 

Figures 34 and 35 show the mass and mixture ratio contours of triplet No. 6 at two 
different operating conditions. Each figure set shows the front and back of the 
contour plot to improve the clarity. At the nomina- mixture ratio this element 
shows evidence of the central oxidizer stream ov< rpowering the fuel streams 
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IRIPLET NO 6 NOM O/F 2IN DI5T. 

Oxldlzar ond Fual Pass Mux (Front sldo) 
Mlxtura Ratio o/f .259262 
Hi King EFf, E* .595469 
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TRIPLET NO 6 NOM O/F 2IN OIST. 
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Figure 34. Triplet #6 Mass Flux Distribution 
at Nominal Test Conditions 
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producing a prominent central oxidizer., mass flux peak, and a relatively low 
mixing efficiency (.595). Reducing the mixture ratio by 20% results in better 
matching of the oxidizer and fuel mass flux distribution, as is evident in the 
improved mixing efficiency (.749). In the review of a new element design, these 
results would suggest a design change in the selected orifice diameters to opti- 
mize element performance at the design nominal mixture ratio. 

Triplet No. 8 . The third and last element to be tested with the liquid/liquid 
mixing system was triplet No. 8. Detail configuration is outlined in Table 7. 
The basic format is a triplet configuration geometrically similar to triplets 
No. 1 and 6. This unit, however, is what is generally referred to as a "reversed 
triplet". This description is applied to impinging elements where the two outer 
orifices are oxidizer streams impinging on the central fuel stream. This element 
is sized for main engine operating conditions with liquid oxygen and liquid RP1 
propellant. The two oxidizer orifices are .165 cm diameter and the fuel orifice 
is .127 cm. The included angle is the traditional 30 degrees each side, with the 
impingement point at .64 cm from the injector face. 

The "reverse" triplet, with nearly the same diameter for both oxidizer and fuel, 
demonstrated a high level of mixing efficiency, and an apparent low sensitivity to 
mixture ratio variations. A "reverse" triplet of similar configuration that was 
tested in an earlier IR&D test series (Ref. 5) indicated somewhat similar charac- 
teristics. The plots of mass distribution also reflect the good mixing efficiency 
in that the shape of the oxidizer and fuel mass distributions are very similar. 

Figure 36 shows the mass and mixture ratio distribution plot for triplet No. 8. 
The oxidizer mass flux (solid line) and the fuel mass flux (dashed line) have sim- 
ilar shapes reflecting the high value of mixing efficiency (93.6%). 

Liquid/Liquid Mixing Analysis of Results 

The liquid/liquid test data have been plotted versus the most commonly accepted 
injector parameters to determine if significant correlation appears to exist. The 
most widely accepted parameter for liquid/liquid triplet elements is the 
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TRIPLET NO. 8, NON O/F 2-IN DIST. 

Ox Idlzar and Fual Ho»» Flux (front »lda) 
Mlxt*a“a Ratio o/F 3.41978 
Mixing Eff/ E» .935961 
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TRIPLET NO. 8. NON O/F/ 2-IN UISJ. 
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Nixing Eff/ Ea .935961 
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Figure 36. Triple! //8 Mass Flux Dist ibutton at 
Nominal Test Conditions 
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Elverum-Morey parameter based on a diameter ratio modified momentum ratio, 
parameter was generated by Elverum and Morey at JPL. Figure 37 shows the liquid/ 
liquid mixing data plotted against this parameter. The purported optimum value 
for a triplet is .66, but this cannot be confirmed with the data obtained to date. 

Also included in this series are plots of mixing efficiency as a function of mix- 
ture ratio, momentum ratio, velocity head ratio, and penetration parameter. None 
of these other parameters is purported to have a best value, with the exception of 
the penetration parameter. The penetration parameter is intended for use with two 
liquid streams impinging on a central gas stream, and a penetration parameter of 
0.5 represents the optimum liquid penetration halfway into the central gas 
stream. The penetration parameter is not generally considered applicable for 
liquid-liquid impinging triplets. 

The data plots for these parameters are displayed as Fig. 37 through 41. Although 
the data is very sparse, there is no obvious correlation between any of these 
parameters and mixing efficiency that applies to all of the three triplet elements. 
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Figure 37. Liquid/Liquid Triplets - Mixing Efficiency 
vs Elverum Morey Factor 
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LIO - LIO TRIPLETS - MIXING EFF. VS SIM. MIXTURE RATIO 
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Figure 39. Liquid/Liquid Triplets - Mixing Efficiency 
vs Momentum Ratio 




RI/RD83-170 


123 


Figure 40. Liquid/Liquid Triplets - Mixing Efficiency 
vs Velocity Head Ratio 
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Figure 41. Liquid/Liquid Triplets - Mixing Efficiency 
vs Penetration Factor 



GAS-LIQUID NIXING TESTS 


Test Results 


Three impinging elements and three coax (or concentric) elements were scheduled 
for the cold flow mixing tests. The impinging injectors (a pentad and two trip- 
lets) appeared to provide reasonable results. The coax element test activity was 
postponted, when relatively extensive testing of coax element 5 indicated serious 
problems. Operation at the simulated mixture ratio conditions indicated very lit- 
tle atomization or mixing of the liquid stream by the gas flow. This character- 
istic indicated problems in the assignment of equivalent flowrates for the test 
conditions. The test facility is not presently capable of matching the gas den- 
sity of the hot-fire case, and is felt that matching either velocity ratio or 
momentum ratio was not truely equivalent for the coax element case. Some effort 
was expended in trying to establish a coax parameter, but indications were that 
the simulated operating points were not satisfactory. The coax element portion of 
this testing has been postponed until the high-pressure facility is available. 

The three impinging elements were tested for a total of 12 useable tests (some 
additional tests were run to establish procedures) . The parameters explored 
included mixture ratio simulations, two different sample distances, and two cham- 
ber pressures (changed gas density). Table 9 shows a detailed summary of element 
tests and test results, and the following section discusses the results obtained 
with each element. 

Triplet No. 2 . Triplet number 2 was conceived as a preburner injector for an 
engine using liquid oxygen and gaseous methane. The design mixture ratio was sig- 
nificantly fuel rich at .49, to provide combustion temperature suitable for tur- 
bine inlet conditions. This was a relatively conventional triplet configuration 
with a central liquid oxygen stream and two gaseous fuel streams (dimensional 
details are shown in Table 7). 

The mixing efficiency of this element in cold flow testing was poor, with the dis- 
tribution plots suggesting that the two fuel streams were overpowering the central 
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liquid stream at all operating conditions. The mixing efficiency in the mid-60's 
does not necessarily reflect low C-star combustion efficiency, but significant 
macroscopic temperature striations would be expected in the combustion process. 
The mixing efficiency did not significantly change with mixture ratio, providing 
no significant clues for establishing parameters for element improvement. One 
test was performed at higher chamber pressure, producing only a slight change in 
mixing efficiency. 

Pentad Wo. 3 . Element No. 3 was a gas/liquid pentad (four-on-one) element 
designed for use in a liquid oxygen-gaseous methane preburner combustor, at the 
same basic operating conditions as the triplet element No. 2. Nominal mixture 
ratio was again .49 for fuel-rich turbine drive gases. A central oxidizer (liquid 
oxygen) stream is impinged by four gaseous fuel streams (dimensional detail is 
outlined in Table 7). 

This element was tested over a mixture ratio range at the baseline test conditions 
of . 025M distance from face to sample plane, and 173 kPa gage (25 psig) sample 
chamber pressure. Two additional tests were run to evaluate a more distant sample 
plane (.05M) and a higher chamber prssure (245 kPa gage). In general this element 
showed high levels of mixing efficiency, with the exception of the data taken at 
greater collection distance. The data taken at the greater distance seems to 
indicate greater over penetration of the gas streams further downstream in the 
chamber. This distance is greater than the anticipated flame front distance, and 
the data at the 2.5 cm plane is judged to be more representative, and was used for 
the majority of gas liquid testing with the impinging elements. 


Triplet No. 4 . Triplet No. 4 also was a preburner element design, differing pri 
marily from element No. 2 by its physically smaller size. This element had been 
designed for a new fabrication process (electrical disposition) which permits 
low-cost replication of very small elements. This also was a relatively conven- 
tional triplet with a central oxidizer stream and two gaseous fuel streams. 
(Dimensional details ar given in Table 7). 
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OF POOR QUALITY 

This element indicated significantly better mixing efficiency than triplet No. 2 
when tested over equivalent ranges. It is difficult to assess whether this higher 
performance should be attributed to its smaller scale, or to the fact that the 
diameter ratio results in higher oxidizer velocity at a given mixture ratio. The 
higher mixing efficiency quantitative data also is supported by the qualitative 
appearance of the mass distribution plot (Fig. 42), which shows the close agree- 
ment between fuel and oxidizer simulant mass flux profiles. 

Coaxial Element Tests . Tentative testing with the coaxial (or concentric) ele- 
ments was very discouraging, to the extent that these elements are being delayed 
to the second phase of testing in the higher pressure facility. The problems 
appeared to stem from the inability of the test equipment to simulate the high 
density of the gaseous fuel (methane) at high operating pressure. Both structural 
and supply system limitations of the gas liquid mixing facility made the simula- 
tion of gas density unattainable. Chamber pressure with this test position is 
limited to about 690 kPa, as compared to operating conditions in the 20.5 MPa 
range for the hot-fire design point. 



FUEL: 

OXIDIZER 

CHAMBER PRESSURE: 173 kPa GAGE (25 psia) 
COLLECTION DISTANCE . .025 M 


Figure 42. Triplet No. 4 NOM MOM R. 
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Attempting to simulate mixing conditions by velocity relationships or momentum 
relationships did not appear to provide the expected dispersal of the central 
liquid stream. The test results appeared to indicate a relatively solid central 
stream of liquid (Fig. 42), which overloaded the collection probe. Operating con- 
ditions were shifted in an effort to establish a rational parameter for simulation 
of the expected hot-fire case, but the results were inconclusive. Increasing the 
gaseous flowrate improved the liquid dispersion and atomization, but the para- 

metric relationship between these test points and the hot-fire conditions were 
tenuous at best. 

Gas-Liquid Mixing Analysis of Results 

The gas liquid mixing data from the impinging elements has been plotted against 
the more common injection parameters, and the results are shown in Fig. 44 through 
48. There is very little information available correlating the gas on liquid 
impinging elements with any flow parameters. The correlating parameters in 
general use are primarily either liquid/liquid parameters or, as in the case of 

the penetration parameter, designed for use with the liquid streams impinging on a 
central gas flow. None of the plots indicate any solid trends with any of the 
selected parameters. The mixture ratio plot (Fig. 44) verifies that these ele 
ments were designed for lower mixture ratios in the preburner or gas generator 
type operating range. The pentad, and smaller triplet were measurably above the 

regular triplet in mixing performance, but no pattern was evident in the data. 

The penetration parameter (nominally devised to describe the penetration of liquid 
streams into a central gaseous stream), Fig. 45, also failed to provide any useful 
trends when applied to the reverse case of gas streams penetrating liquid. The 
Elverum-Morey parameter (Fig. 46) also was disappointing, with each group of data 
seeming to start its own curve, with no interaction with the curves from the other 
elements . 

The other parameters (momentum ratio and velocity head ratio. Fig 47 and 48), also 
show little tendency to provide a single orderly curve. The most likely con- 
clusion from this is that the elements are too varied in geometry and design to 
produce a single curve. Serious parametric studies would require working with a 
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Figure 44. Gas/Liquid Impeller Elements - Mixing Efficiency 
vs Simulated Mixture Ratio 
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Figure 45. Gas/Liquid Impeller Elements - Mixing Efficiency 
vs Penetration Factor 
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Figure 46. Gas/Liquid Impeller Elements - Mixing Efficiency 
vs Elverum Morey Factor 
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Figure 47. Gas/Liquid Impeller Elements - Mixing Efficiency 
vs Momentum Ratio 
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LIQ IMP ELEMENT? - MIXING EFF. VS VELOCITY HERD RRTIO 



single family of elements, varying a minimum number of characteristics between 
elements. Minor manufacturing discrepancies appear to be capable of producing 
significant shifts in mixing characteristics. 


The disappointing results of the coaxial (or concentric) elements resulted in 

reduced confidence for all the gas-liquid tests because of the inability to match 
reactant density values with the simulant fluids. The gas liquid mixing tests 

were hampered by facility limitations in simulating the gas injection density. 
Most of the elements selected for this test series were designed to operate with 
gaseous hydrocarbons at high chamber pressures. These gaseous hydrocarbons are 

many times denser than the hydrogen or hot hydrogen/oxygen combustion products, 
for which most previous data was obtained. The gas-liquid mixing chamber used in 
these tests was structurally limited to about 690 kPa (100 psig), and in most of 
these tests, the supply systems limited test operation below this value. With 
impinging element injectors, the most widely accepted correlation factors are 

based on the ratio of the stream momentums . With the diameter ratio fixed by com 
mon hardware matching momentums ratio is the most accepted method of simulation 
with fluid densities differing from the hot-fire conditions. In the liquid/liquid 
mixing tests, this was a minor correction, since the simulant density ratio 
quite close to the reactant density ratio. For gas/liquid test this density ratio 

is significantly different. 

The density of injected gaseous methane is on the order of 160-190 kg/M 3 (10 to 
12 pounds per cubic foot) when the methane is employed as an engine coolant prior 
to injection. The maximum nitrogen gas density obtainable in the cold flow test 
apparatus is on the order of 9 kg/M 3 (.57 pounds per cubic foot), or a factor of 
20 lower. Oxidizer to fuel density ratio for hot fire is typically 6 to one, 
while the lowest value obtainable in the available test hardware is 110 to one 
(the higher pressure test chamber now being pre pared will be a significant 

improvement) . 


Matching the design point momentum ratio to simulate the reactive test conditions 
does not automatically match the other parameters. Those parameters that vary 
with the density ratio of the flowing fluids will end up with different values 
than for the hot-fire case. The various plots of mixing efficiency as a function 
of various parameters (Fig. 44 through 48) arc plots of the cold flow parameters, 
using cold flow simulant mixture ratio, density ratio, etc. 
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MIXING TEST CONCLUSIONS 


Both liquid-liquid, and gas-liquid, cold flow mixing tests were performed and the 
resultant data reduced for analysis. Three different triplet elements were tested 
in the liquid-liquid program. Two triplet elements, a pentad and a coaxial ele 
ment, were tested in the gas-liquid portion of the evaluation. Two additional 
coaxial elements were scheduled for test, but the limitations of the test equip- 
ment has resulted in postponement of these tests, until a higher pressure facil- 
ity, now being prepared, becomes available. 

The liquid-liquid elements tested consisted of two elements sized for fuel- rich 
operation in gas generator or staged combustion preburner applications, and a 
’’reverse” triplet sized for main chamber operating conditions. The two elements 
designed for low-mixture ratio exhibited depressed mixing efficiency due to the 
unbalance in flow and diameter ratio, while the reverse triplet indicated good 
mixing. Plotting the data from these three elements did not appear to support any 
continuous curve of mixing efficiency versus any of the accepted mixing parameters. 


The gas-liquid mixing tests encompassed two triplets, and a pentad, impinging 
designs, and one coaxial injector. Three coaxial elements had been fabricated, 
but the results of the first coaxial element to be tested indicated the critical, 
need for higher pressure testing. Therefore, the coax part of the program has 
been postponed until the high pressure facility is ready. 

All three gas/liquid impinging pattern elements were sized for low-mixture pre- 
burner or gas generator operation with gaseous fuel and liquid oxygen. One of the 
triplets is a micro-orifice injector, which is characterized by a high-element 
density (elements per square inch) and small orifice sizes. 

One of the goals of this program is to establish, or reinforce, injector correla- 
tions between mixing characteristics and injector design parameters. The primary 
benefit from this type of result would be to provide design guides for the selec- 
tion and sizing of injector elements for new applications, and the ability to 
predict performance parameters of new designs prior to costly fabrication and test. 


RI/RD83-170 

137 



The results of the testing to date have not been particularly promising from this 
standpoint. These test results do not provide any high degree of correlation with 
any of the existing mixing parameters. Several factors undoubtedly led to this 
condition and point out the difficulty in providing any sort of quantitative 
predictive parameter for a wide range of injection phenomena. Two of the elements 
tested had a significant unbalance in injection orifice diameters, as a result of 
operation at the extremes of mixture ratio associated with preburner and gas gen- 
erator operating conditions. There also were some evidence of manufacturing or 
manifolding problems upsetting the symmetry of these elements. 

To establish a proven injector mixing correlating design parameter requires a sys- 
tematic program with duplicate configurations and redundant testing to overcome 
statistical spread in test results. In this manner, the random variation in 
geometry (e.g., misimpingement , smoothness of holes, hole entrance conditions) and 
flow conditions, which can be minimized but never eliminated, will be "averaged" 
out . 

Mixing tests establish the degree of "mechanical" mixing (i.e., the mixing that 
occurs in the absence of combustion and combustion-induced gas motion). As such, 
they provide only a rough indication of the actual degree of mixing in an engine. 
Nevertheless, optimization of the mechanical mixing is a reasonable prerequisite 
for any injector design. Thus, mixing correlating parameters that can be utilized 
as design tools are an important technical need. The literature survey performed 
as a part of this effort demonstrates the lack of information regarding the mixing 
charactristics of, and the lack of verified mixing correlating parameters for, 
triplet, pentad, and coaxial elements. Thus, the need for a research program such 
as this is apparent. 

The cold flow mixing testing has experienced a number of problems. Some of these 
are still being resolved and some serve as "lessons learned" and issues to be 
resolved in future testing. These problems and issues include (1) gas/liquid test 
collection efficiency deviations, (2) element imperfection effects (differences in 
mixing characteristics produced by "identically" constructed elements), (3) col- 
lection distance effects, (4) the need for high pressure testing, and (5) the 
lack, to date, of verification of the existence of, and the optimum values of 
mixing correlating parameters (especially for coaxial elements). 
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The cold flow test hardware used in low-pressure mixing tests is presented here. 
Drawing 7R001670 depicts the gas/liquid and liquid/liquid element test specimens, 
as well as the necessary manifolds and mounting hardware. 
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APPENDIX B 


This appendix describes the manner in which the raw test data is converted to the 
presented test results. The means by which mixing is measured (primarily mixing 
efficiency) are defined, and the manner and software by which these means are com- 
puted are described. The computer programs and examples of their input and output 
are presented. Data reduction methods for liquid/liquid and gas/liquid mixing 

tests are reported separately. 

LIQUID-LIQUID DATA REDUCTION 

The liquid-liquid mixing test data is hand recorded on a special form that aids in 
entering the information into microcomputer. The data for each sample tube xs 
recorded in cubic centimeters as the total sample and the "Trie" sample level (the 
heavier trichloroethane settles to the bottom of the sample tube, and the coloring 
agent helps provide the demarcation line). A data file is typed into the micro- 
computer (Xerox 820) and saved on "floppy disk" for each of the test data sheets. 
This data file includes a text line describing and identifying the test, and 
several lines of overall test information with such things as pressures, flow- 
rates, test times, element injection orifice diameters etc. The detail sample 
tube data is entered in a tabular format according to sample location as copied 
from the test data sheet. A sample of the test data sheet is presented as Fig. 

B-l. 

The data reduction program was converted and updated from programs used in the 
past, and is configured to run in Micro soft BASIC on the Xerox 820 microcomputer 
(listing and sample data file are shown in Fig. B-2). This program converts the 
volumetric sample sizes in the input to equivalent propellant local flowrates in 
"stream tubes." This conversion is based on simulant and propellant densities. 
If the density ratio of the simulants is fairly close to the reactants, this cor- 
rection factor is relatively small. The simulants, trichloroethane at 1320 
Kg/M 3 , and water at 1000 Kg/M 3 , have a density ratio of 1.32 which is not too 
far from the "L0X" (at 1137 Kg/M 3 ) and RPI (at 800 Kg/M 3 ) ratio of 1.42. The 
liquid oxygen/liquid propane density ratio of 2.3 is further from the simulated 


value . 


RI/RD83-170 

B-l 



The most important factor in cold flow mixing data reduction is, of course, the 
determination of mixing effectiveness. The parameter that is generally accepted 
for this determination is the mixing efficiency, frequently referred to as E sub 
M”, the use of which was pioneered for rocket injectors by Jack Rupe at JPL. This 
mixing efficiency is computed by a mass weighted summation of all the samples, 
decrementing each sample as a function of its local oxidizer fraction related to 
the overall oxidizer fraction of the total sample. 


The oxidizer fraction (oxidizer sample weight over (he total sample weight) is 
used in place of the more usually referenced mixture ratio (oxidizer weight flow 
over fuel weight flow) because it is a more rational mathematic relationship, 
ranging from zero to one rather than zero to infinity. The computer program com- 
putes mixing efficiency based on the local oxidizer traction and the local mass 
collected relative to the overall total collected oxidizer fraction and total mass 
collected . 


The nominal form for computing mixing efficiency is: 



L 


where 

0 



M 


sa 


0 


sb 


0 


sa 


Mixing efficiency. Maximum = 1 00V 

Overall oxidizer fraction expressed as total oxidizer weight 
over total sample weight 

Mass fraction (total weight of local sample over total 
weight of all samples) subscript sb defines those samples 
with oxidizer fractions below the overall oxidizer fraction 
value 

Mass fractions of samples with local ox fractions above the 
overall value 

The sample oxidizer fraction value for those samples below 
the overall value 

The sample oxidizer fraction value for those samples above 
the overall value 
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The mixing efficiency of a set of samples, all at the same oxidizer fraction 
value, would be 100%, regardless of how the sample mass is distributed (so long as 
the mixture ratio is uniform). The opposite condition, where all the oxidizer (or 
oxidizer simulant) is in one set of samples and all the fuel in another set, will 
result in mixing efficiency of zero. Mixing efficiency values are much more sen- 
sitive than engine combustion efficiency, and it is not unusual that a mixing 
efficiency in the mid-80% range may indicate combustion C-star efficiencies that 
are much higher. 


Frequently in evaluating mixing data, it is useful to utilize a simple stream tube 
analysis to predict mixing effects on measured combustion efficiency. The most 
common approach to this task is to assume that each sample represents a "tube" of 
the combustion process with little or no cross mixing between sample "stream 
tubes". Full-efficiency combustion is assumed for the mixed products within each 
"tube”. The resulting C-star is multiplied by the total mass flow in the sample. 
These values are added, and the product is divided by the total mass flow of all 
the samples. The result is an approximation of the C-star resulting from the 
mixed flow. The form of this computation is: 


C* . 


mix 



^ theo(l) + ^(2) ^ theo(2) 


M 


tot 



C* 



where 


c * 

M, 


mix 


(1) 
C *theo(l ) 
M tot 


= Mixing Limited C-Star 
= Mass flow in respective samples 
= Theo C-star at local sample mixture ratio 
= Total collected sample mass flow 


The mixture limited C-star computations are generally very optimistic for 
operating conditions at non-stoichiometric mixture ratios. In these operating 
ranges, the low— mixture ratio "stream tubes’ are balanced (and sometimes 
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overbalanced) by the high-mixture ratio samples. Two of the liquid/liquid ele- 

ments tested in this program are designed for very low-mixture ratio operation in 
gas generator or preburner applications. The theoretical C-star curves in these 
low-mixture ratio regions cause very strange results for this type of computation, 
frequently predicting over 100% mixing limited C-star efficiency, even with poorly 
mixed elements. For this reason, the mixing limited C-star computations are not 
being utilized in this program. 

The collection efficiency is the quantity of fluids collected divided by the quan- 
tity of flow through the injector during the test. Values far from one indicate 
some error or problem with the test procedure. For liquid/liquid tests, col- 

lection efficiency can be inferred from the test data in the following manner. 

Collection efficiency is a method of cross checking the test sample against the 
theoretical sample element flow during the test period. With the liquid/liquid 
flow facility, any sample that escapes the collection grid can be inferred from 
the shape of the mass distribution plots. If the mass distribution plots show 
significant mass collected near the edge of the collection grid, this is strong 
evidence that measureable sample spray fell outside the sample grid. The quanti 
tative measure of collection efficiency will assist in the decision of whether a 
rerun of the test with a repositioned element, is warranted. If the spray appears 
to have exceeded the collection grid, collection efficiency can be computed in the 
data reduction program by summing the collected samples and comparing to the 
computed value of total flow through the model element during the sample time. 
Generally this computed total flow is based on the calibration flow resistance of 
the element model and inlet pressure. For the liquid/liquid tests of this 
program, the collection efficiencies were very near to one. 
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DATE: 11-9-82 

RUN #: +20% MR-MOM RAT 

INJECTOR#: TRIPLET 1 



/• til 

QUALITY 


PSID (WATER): 17 + 1 = 18 
PSID (1 ,1,1,-TRIC) : 36 + 3 = 39 
TIME DURATION: 15 SEC 


COLLECTION DISTANCE: 2" 



0 0 .25 .25 0 0 

0 1 3.5 3 2.5 .5 



2 3 4 5 6 7 8 9 10 11 12 


Figure B-l. Liquid/Liquid Mixing Sample Data Sheet 
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1 REM PROGRAM CI.DFLTC I. IQ/I Q COI I) FI 0/ MIXING DATA REDUCT 

2 REM UPDATED VERSION FEB 28 83 

3 REM VERSION FOR TRIPLETS - NORMAI AND REVERSED 

4 REM 3D DATA PLOT OR PI OT FILE OPTION 

5 REM PROGRAM CLDFLTC LI 

7 REM WITH ROTATED DISPLAY AND PEN 1 FOR CX PEN 2 FOR FU 

8 REM READ TEXT 

9 REM READ INJ NO, RUN NO,TRIC DELP,I120 DEI P, COI I TIME 

10 DATA "TRIPLET NO 1 -+20%MXR - 2 IN" 

11 DATA 1,1,39,13,15 

12 DATA .058, .71 ,2, .0465, .814,1 

110 DATA 0,0, 0,0, .5, .5, .25, 0,0, 0,0, 0,0 

111 DATA 0 , 0 , O,0, 4, 5, 3, 1,0, 0,0, 0,0 

120 DATA 0, 0,0,0, .25, .5, .5,0,0, 0,0,0, 0 

121 DATA P,0, 0,0, 4, 6,4,1. 5, 0,0, 0,0,0 

130 DATA 0, 0,0,0, .25, 1 , .5,0,0, 0,0, C,0 

131 DATA 0,0, 0,0, 3, 7, 6, 1.5, .5, 0,0, 0,0 

140 DATA 0,0, 0,0, .25, 1.5, 1,0, 0,0, 0,0,0 

141 DATA 0,0, 0,0, 3, 9, 8, 2, .5, 0,0, 0,0 

150 DATA 0,0, 0,0, .25, 2, 2, .25, 0,0, 0,0,0 

151 DATA 0,0,0,0,2.5,10,11,3,1,0,0,0,0 

160 DATA 0,0, 0,0, .25, 2, 2. 5, .25, 0,0, 0,0,0 

161 DATA 0,0, P,0, 2, 11, 16, 5,1, .25, 0,0,0 

170 DATA O, 0,0,0, . 25, 2,4.5, 1 , .25,0, P, 0, 0 

171 DATA 0,0, 0,0, 2, 11, 24, 7, 1.5, .25, 0,0,0 
180 DATA 0, 0,0,0, .25, 2, 6, 2. 5, .25, 0,0,0," 

1 3 1 DATA 0,0, 0,0, 2, 11, 31, 11, 2, .25, 0,0,0 
190 DATA 0, 0,0,0, .2 5, 2,7. 6, 3, .5, 0,0, 0, 

■ r 


233 

DATA 

9,9,3 

,0, .25,2.5,8,5, .5, 0,0,n 

2 31 

DATA 

0 ,0,9 

, 9, 2, 19, 34, 29,4,. 5, i:, 9,; 

213 

DATA 

3,0,0 

,9,. 25,/ *5, 3, 5, 1,0,1 ,3,1 

211 

DATA 

0,3,9 

,0,2, 0,23, 20, 5, .5, 0,0,0 

223 

DATA 

3,0, 0 

,0, .25,2.5,8,6, 1 .5, .25, 

221 

DATA 

0,0, 3 

, 0 , 1 . 5 , 7 , 2 1 , 1 ? , 5 , 1 , 3 , 9 , 1 

2 33 

DATA 

8,0,0 

,0, .25,2,7.5,6, 1,0, 0,0,1 

231 

DATA 

0 ,0,0 

,0,1,5,19,14,4,1 ,P,F,P 

243 

DATA 

9 ,3,3 

r 0 t 9 t 1 § 5 t 4 r • 5 , / >1 1 0 t 0 

241 

DATA 

3 , i: , 9 

r 9 , .5,4,1 3, 19,2.5, • 5 , 3 # 

2 5 3 
251 

DATA 

DATA 

9,3,3 

, 3 , 9 , . 2 5 , 2 . 5 , 2 , . 2 5 , 9 , 3 , 

'3 ,3,9 

’, O ’, 2 . 5 , 8 , 7 , 1 . 5 , . 5 , 0 , 0 

263 

DATA 

3,9,9 

,0,0, 8,1,1, .25, 0,8, 0,0 

261 

DATA 

3 , 9 , 3 

,8,0, 1,5, 6, 2. 5, .5, 0,0,0 

27 3 

DATA 

3 ,9,9 

, 0 , 0 , 0 , . 5 , . 5 , . 2 5 , 0 , 0 , 0 , 

271 

DATA 

9,3,9 

, 9 f ( , .5,3,5, j,l,..,i , 

283 

DATA 

0,9, <9 


283 

293 

DATA 

r» r* r» 

, 9 , 0 , 9 , 1 , 3 . 5 , 3 , 2 . 5 , . 5 , 9 

DATA 

d/i. 

,0,0, 0,0, .75, .25,0,0,0, 

291 

DATA 

3,0,9 

,0,0,1 ,.5, 2, 2.5, 2. 5, 1,0 

333 

DATA 

9,3, 0 

,0,0, 0,0, .25,0,8,0,0,0 

301 

DATA 

n ^ :« 

,8, 0,0, 0,1, 1,1. 5,1, .25, 

19 r 

9 REM 

COI D 

FIOW -BASIC VERSION OF 

1923 n EM 

WDOT 

IN IB MASS/SEC DENSITY 

1030 DIM 

TRIC ( 

20, 1 3) 

1343 DIM 

C’X S M P ( 2 0 , 1 3) , TOT (20,13) , H70 ( 


, OXERS M ( 2:’ , 1 3) 


Figure B-2. Liquid/Liquid Mixing Data Reduction Program 
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QUALITY 


1050 DIM FUSMP(20,13) , TOTP (20,13) , XMR (20,13) , DEC (20,13) 

1060 DIM FUM (20,13) , OXUM (20,13) 

1070 DIM OFRTHE(10) ,CSTIIE(10) 

1080 DIM TMTHE(IO) 

1090 DIM SPIMPS(IO) ,SPIMPF(10) ,CSTHES(10) 

1108 DIM TOTS (20,13) 

1102 PRINT" INPUT DATE AND TIME - ANY FORMAT" 

1103 INPUT RTS 

1110 PRINT" INPUT CODE FOR OUTPUT 0 SHORT, 1 DEI INP, 2 FUM. OUTPUT" 

1120 INPUT I PRT 

1122 IPRINT"DATA REDUCTION RUN OF " ; RT$ 

1130 REM ZERO OUT VARIABI.ES 
1140 CSSUM-0 
1150 DECSUM-0 
1160 TOTPSM-0 
1170 TOTOSM-O 
1100 TOTFSM-0 
l l'.'0 DENTRC-82.59 
17. ' . 'K i'= '« 7 . 7 7 
./! • j = 

dkmox-71 

1230 SAMPNO’260 

1 1240 FUCOR-. 0022* ( (DENFU/DENWAT) * . 5) 

1250 OXCGR*. 0022* ( ( DENTRC/DENWAT) * ( ( DKNOX/DENTRC) “ . 5) ) 

1260 REM 

1278 REM READ DATA 
1200 READ texts 
1290 PRINT TEXTS 
13 00 I PRINT 
1310 I PRINT TEXTS 
; • r Pt» i 

i ’ ! 1 i vj , ite::t , p i ti< ic , pi 1120 

1 340 REM INJ NO, TEST NO, INI E T PR TR1C, 1120 

1 750 I PR 1 NT" I N J " ; INJ;" TEST NO. " ; I TEST; " INI ET PR TUIC ";PITRIC,*" I12U ";PIM20 
1360 READ TIME 

1370 READ DWAT, C DWAT , WATN , DTR IC , CDTU IC , TR ICN 

13P0 REM DIA, CD, AND NO. OF WATER AND TFIC INJ ORIFICES 

1 390 [PRINT" OX DIA " ; DTR I C ; TR ICN ; " PI ACfc!> FU DIA " ; DWAT ; W ATN ; “ PEACES" 

1400 [PRINT 

1410 AREAFU* .7854 *V.ATN* ( DWAT* 2 ) 

1420 AREAOX* .7054 *TR ICN* (DTRIC‘2> 

1430 FOR KROW-1 TO 20 
1448 FOR KCOI»l TO 13 
1450 READ TRIC (KROW, KCOI ) 

1460 IF TRIC (KROW, KCOI )<. 000.7001 THEM TUIC ( KROW, KCOI )». OO0UOO1 
1470 NEXT 

1430 FOR KCOI *1 TO 13 
1490 READ TOT (KROW, KCOI.) 

1500 IF TOT (KROW, KCOI X. 0001 THEN TOT ( KRCW , KCOI. ) » . Oil 8 1 

1510 NEXT 

1520 NEXT 

1530 REM 

1540 REM 

1550 IF I PRT< 2 THEM 1700 
1560 REM DIAGNOSTIC PRINT OUT 


Figure B-2. (Continued) 
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1570 L PRINT" INPUT DATA CC - TOT / TRIC" OF FOOfi QU-^^i 3 

1580 L PRINT 

1590 FOR KROW-1 TO 20 

1600 FOR KCOl-1 TO 13 

1610 L PRINT USING”##.# " ? TOT (KROW , KCOI ) ; 

1620 NEXT 

1630 LPRINT 

1640 FOR KCOL-1 TO 13 

1650 LPRINT USING"##.# " ; TRIC (KROW , KCOL ) ; 

1660 NEXT 
1670 LPRINT 
1680 LPRINT 
1690 NEXT 

1700 REM BYPASS W/G PRINT 
1710 REM 
1720 REM 

1730 REM INITIAriZE MATRIX SIZE 

1740 NROWS-1 

1750 NROWE-28 

1760 NCOIS-1 

177(1 NCOLE* 1 3 

1780 TOTPSM-C 

1790 TOTTRC- 3 

1800 TOTS IM-C 

1310 TOT0SM-0 

1320 TOTFSM-U 

1830 TOTWAT-C 

1340 NMSAMP-0 

1850 IF I PRT-2 THEN [PRINT 

I860 TF I PRT- 2 THEN [.PRINT" MIXTURE RATIO" 

1970 IF T PRT-2 THEN [PRINT 
133 0 FOR K R OW * N R OW S TO N R OR <E 
1890 FOR KCOI-NCOIS TO NCOl E 

1<>00 IF TOT (KROW , KCOI ) < TRIG { KUO a , KC 01 ) THEN TOT ( KR'"\ , KC Of ) - . 0 01 +TR TC ( K ROW , r ' 

1 9 1 0 1120 ( KROW , KCOI ) -TOT ( KROW, KCOI ) - TR I C (KUO* , KCOI ) 

1920 REM 

19*8 REM TEST FOR ZERO AMOUNTS 

1940 IF TOT (KROW, KCOr ) > . 35 THEN N M 3 AM p - N % l SAM I 5 + ] 

1 960 REM 
I960 REM 

1970 IF H20 (KROW, KCOI ) *0 THEN 2880 

1 "8 3 XMR ( KROW , KCOL ) -TRIC { KROW , KCOI ) /!! 20 ( KROW , KCOI ) 

199 8 XMR (KROW, KCOI ) -XMR (KROW , KC 01 ) *OXCOR/FUCCR 
28C8 REM 

2810 IF IPHT-2 THEN IPRINT USING"##. H !| H ; XMR ( KROW , KCC I ) ; 

2 020 OXSMP (KROW, KCOL) -OXCOR *TR TC ( KROW , KCOL ) 

2 0 38 FUSMP (KROW , KCOI ) -FUCORM120 ( KROW , KCOI ) 

2 04 8 TOTP ( KROW , KCOI ) -OXSMP { KROW , KCOI ) +FUSMP ( KROW , KCOI ' 

2 0 50 TOTS (KROW , KCOI ) -TRIC ( KROW , KCOI ) + !120 (KROW , KCOI ) 

2060 FACTOR- I 

2070 TOTP {KROW, KCOI ) -TOTP ( KROW , KCOI ) +FACTOR 
209 8 TOTS (KROW, KCOL) -TOTS (KROW, KCOI ) * FACTOR 
2090 TOTOSM-TOTOSM+OXSMP ( KROW , KCOI ) 

2100 TOTPSM-TOTPSM+TOTP (KROW, KCOI ) 

2110 TOTTRC -TOTTRC +TR IC (KROW , KCOI ) 


Figure B-2. (Continued' 
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2120 TOTSIM-TOTSIM+TOTS (KROW,KCOL) 

2130 TOTFSM-TOTFSM+FUSMP (KROW, KCOI ) 

2140 T0TWAT-T0TWAT+H20 (KROW, KCOL) 

2150 OXFRSM (KROW, KCOL) -OXSMP (KROW, KCOL ) /TOTP ( KROW , KCOL ) 

2160 NEXT 

2170 IF IPRT-2 THEN LPRINT" " 

2180 NEXT 
2190 LPRINT 
2200 REM 

2210 REM 101 CONTINUE 
2220 REM OVERALL SUMMATIONS 
2230 XMROA-TOTOSM/TOTFSM 
2240 OXFOA»TOTOSM/TOTPSM 

2250 REM COPUTE FLOWRATES FROM INPUT PR,DIA,CD 

2260 WWAT»CDKAT*SQR ( DENWAT* PIII20* ( WATN * 2 ) * ( DWAT* 4 ) /3 . 625 ) 

2270 WTRC»CDTRIC*SQR(DENTRC*PITRIC* (TRICN“ 2 ) * (DTR IC* 4 ) /3 . 625) 

2280 WFUI?J»WWAT*SQR ( DENFU/DENWAT) 

2290 WOXIN»WTRC*SQR (DENOX/DENTRC) 

2300 WWATC»3. 5 3 1 E-0 5 *TOTVvAT* DENWAT/T IME 
231 0 WTRCC* 3 . 5 3 1 E-0 5*TPTTRC*DENTRC/TIME 
2320 OFRIN-WOXIN/WFUIN 
2330 OXFRIN*WOX IN/ ( WFUI N+WOX IN ) 

2340 OXTI N*TIME*W'OX I N 
2350 FUTIN»TIME*WFUIN 
2 360 OXCI EF-WTRCC/WTRC 
2 370 F'JCf EF-WWATC/WWAT 

2380 [PRINT" COLLECTED OX SIM ";TOTOSM;" FUEI SIM "-T0TFSM 
2390 I PRINT" COMPUTED OX SIM ";0XT1N;" FUEr SIM "-FUTIN 
2400 I PRINT" COMPUTED COI D F!UW O/F " ; (WTRC/WWAT) ' 

2410 LPRINT" COMPUTED SlMUTATED O/F " ; (WOXIN/WFU1N1 
2420 I PRINT "COII EFF OX ";OXCLEF;" FU ";FUCIEF 
2430 I PRINT" COIIECTED TRIC/WAT O/F ” ; (WTRCC/WWATC ) 

24 50 REM INT " COn ' ECTED SIMUI ATt:D O/F " ; XMROA; " OX FRACT ";OXFOA 

2460 REM E SUB EM COMPUTATIONS 
2470 TOTPST-O 
2480 CSSUM»0 
2490 DECS UM» 0 

25011 FOR KROW»NRCW5 TO NR OWE 

2510 PRINT" ROW " ; KROW 

2520 FOR KC 0! =NC OLS TO NCOI E 

2530 IF OXFRSM ( KROW , KCOI ) >» OXFOA THEN DECR = OXFRSM ( KRO’.v , KCOI ) -OXFOA 

2540 IF OXFRSM (KROW, KCOL) < OXFOA THEN DECR*OXFOA-OXFRSM ( KROW KC^I 1 

2550 DECRT»DECR*TOTP (KROW, KCOI ) ' 

2540 DECRT«DECKT/OXFOA 

2570 DECSUM-DEC3UM+DECRT 

25C0 T0TPST*T07P5T+TGTP( KROW, KCOI ) 

2590 REM 

2600 REM MIXING I IM I TED CSTAR 
2610 REM INPUT MIX IIMIT CSTAR CAICS 
2620 REM 
2630 REM 

2440 REM 331 CONTINUE 
2650 NEXT KCOI 
2660 NEXT KROW 


Figure B-2. (Continued) 
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2670 ESUBM* (TOTPST-DECSUM) /TOTPST 
2630 LPRINT" E SUB M EFF " ; ESUBM 
2690 CSTAR-CSSUM/TOTPST 

2700 REM COMPUTE THEO CSTA: : < O/AI I OXRATIO 

2710 REM COMPUTE CSTAR EFF 

2720 REM INJECTOR PARAMETERS 

2730 REM CONVENT IONA I TRI PIETS 

2740 WOUT-WWATC 

2750 WIN-WTRCC 

2760 DENOUT*DENWAT 

2770 DENIN-DENTRC 

2733 AREAOT-AREAFU 

2790 AREAIN-AREAOX 

2800 DIAOUT-DWAT 

2810 DIAIN-DTRIC 

2820 I PRINT 

2333 IF WATN>TR ICN THEN 2950 

2840 REM REVERSE TRI PIETS 

2850 I. PRINT" REVERSE (O-F-O) TRIPIET” 

2860 WOUT-WTRCC 
2870 WIN-WWATC 
28 8 0 DENOUT=*DENTRC 
2890 DENIN-DENWAT 
2900 AREAOT-AREAGX 
2910 AREAIN-AREAFU 
2920 DIACUT-DTRIC 
2930 DIAIN*DWAT 
2940 GOTO 2970 
2950 REM 

2960 t.PR INT" NORMAI (F-G-F) TRIPIET" 

2970 REM 
2983 REM 

2990 V E I R AT* WTRCC * AR E A F U * D E NW AT / ( W'.\ ATC* AH EAOX*DFNTRC ) 

2992 VEI HDF-2.23S* (i*WATC"2) / (DENWAT* (AREAFU‘2) 1 
2994 VEI.HDO-2.236* ( WTRCC* 2 ) / ( DENTRC* (AREAOX'2) ) 

2996 VEI HDR-VEI HDO/VEI liDF 
3000 MOMRAT»WTRCC*VEI RAT/WWATC 
3010 REM El V MOREY 

30 2i’ ei i.vA-: = ( (•■ • . . . .;.)*(( v;;ka i n/a.rea; rj * 1 .76} 

3030 REM PRINT OUT INJ PARAMETERS 
3040 I.PRINT 
3050 I PRINT 

3060 I PRINT"CCI I.ECTED CONDITIONS- INJECTION PARAMETERS 
3070 I. PRINT 

3080 l PR INT" MOMENTUM RATIO ";MOMRAT 
3090 I. PR I NT" VEI RATIO " ; VEI RAT 

3092 I.PRINT" VEI MEAD OX SIDE " ; VEI 1 1 DO; " PS 10" 

3094 I PRINT" VEI HEAD FU SIDE " ; VEI 11DF; " PSID" 

309r', I PR I NT" VEI HEAD RATIO "jVEIHOR 

3100 I PRINT"EI VERUM-MOREY FACTOR " ; E t MRA7 

3113 REM PENETRATION PARAMETER ( 3 0 DC, REE ASSN. ED ANSI) 

3120 PENT*1 . 25*SIN (30/57.3)* (WOUT/WIN) * ( ( DEN I N/DENOUT) * . >) * CM A T N/P IACJT ) 
3130 I PR I N’T" PENETRAT I ON PARAMETER ";PEUT 
314 0 I. PR INT 

3150 REM CALCUI.ATE NORMAI I ZED MASS FI UX FU AN OX 


Figure B-2. (Continued) 
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ii no 
M 7i1 
1 1 s :i 
3190 
3200 
3210 
322 0 
3230 
3240 
3250 
1250 
3270 
3230 
1290 
3100 
3310 
3320 
3330 

134 0 
3351? 

135 0 
3370 
3 3 B 0 
1390 
1 4 W 0 
1410 
3 4 2 0 
14 30 


OF Rv 

FOR KHOW-1 TO 20 
FOR KCOI-1 TO 13 

F UM (KROW, KCOI.) -FUSMP (KROW, KCOI.) *NMS AMP/TOTFSM 
OX UM ( KROW, KCOI. ) ~OXSMP ( KROW , KCOL ) *NMSAMP/TOTOSM 
NEXT KCOI. 

NEXT KROW 

REM WRITE FILE FOR PLOT PROGRAM 

PRINT" IF YOU WANT 0 WRITE PLOT FILE INPUT Y H 

INPUT PL $ 

IF PL$<> " Y" THEN 3440 

PR INT" INPUT DESIRED NAME OF PLOT FILE" 

INPUT PFIL$ 

PF I L$*PF I L $ + " , BAS " 
open" o" , *i,PFir.$ 

LN*10 

PRINTS 1 , T N ; "DATA " ; TEXT 
LN-LN+10 

PR I NTS 1 , LN ; "DATA " ; XMROA 
I N*I N+l 0 

PR I NT #1 , LN ; " DATA ";ESUBM 
FOR KROW*l TO 20 
FOR KCOr-1 TO 13 
I.N»IN + 10 

PR I NT# 1 , I N; "DATA 

WRITES 1 f OX UM (KROW, ICC U I ) , F UM { KROW', KCOI ) 

NEXT 
NEXT 
Cl OSE 



3440 REM PLOT ROUTINE FOR MASS DIST 

1450 PRINT” IF YOU WANT TO PI OT MASS DIST PI tJG IN PI OTTER AND ■ 
3400 PRINT"GET THE PEN AND PAPER READY THEN INPUT Y” 

3470 INPUT PTOT$ 


34R0 IF PLOT$< > M Y M TilEN 4140 
349 0 I PRINT" ; : IOD 0 H A PI U " 

3500 XPO»13flfl 
3510 YPO* 300 
3520 VSCAI E~100 

3530 r PRINT” 300,200 SI 2 PLOT OF OX AND FUEL MASS FI UX 
3 54 0 L PR I NT" 3 50,1 50 SI 2 " ;TEXT$;" _ " 

1550 I PRINT” 3 50,100 S12 MIXTURE RATIO XMROA;" " 

35 CO I PRINT" 350,50 SI 2 MIXING EFF - L SUB M " ;ESUBM; M 
3570 REM PLOT FUEI DIST 
3530 I PRINT" P2 " 


3590 FOR KROW-1 TO 20 

XI PO»XPO-30* ( KROW - 1 ) 

** ! n i - v r j - f-so* (ut-,-1 ) 

152 i \:\ j;cei =1 1 , 1 > 

3 5 30 KCO*KCOI-l 

3540 IF KCO< 1 THEN KCO-1 

555 0 XOFF-XLPO+S0* (KCOI -1 ) 

3 5 5 0 YCFF*YI PG+15* (KCOI -I) 

3570 XPT *INT (XOFF) 

3 510 YPT * I NT ( FUM ( KR OW , KC 01. ) * VS C A I. E + Y 0 F F + . 5 ) 

15 9 0 I PRINT" 10 " 

3 70 0 IF (FUM (KROW , KCOI ) + FUM (KROW, KCO) ) < (OX'JN (KROW, KCOI ) + OXUM ( KROW , KC 0 ) ) 


THEN LPRINT" r 


Figure B-2. (Continued) 
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3710 [.PRINT XPI ; " , " ; YPI. ; " D" 

3720 NEXT 

373 0 I. PRINT" U " 

3740 LPRINT" 10 " 

3750 NEXT 

3760 REM PLOT OX DIST 
3770 I PRINT" P3 " 

3780 FOR KROW-1 TO 20 
7790 XI PO-XPO-30* (KROW-1 ) 

7800 YI PO-YPO+50* (KROW-1 ) 

3010 FOR KCOI »1 TO 13 
3820 KCO-KCOI-1 
3830 IF KCO< 1 THEN KCO-1 
3840 XOFF-XLPO+83* (KCOI.-l ) 

3850 YOFF-YLPO+15* (KCCI.-l) 

3560 XPT -INT(XOFF) 

3370 YPI -INT {OXUM (KROW, KCOI ) *VSCAI. E+YOFF+ . 5 ) 

3380 IPRINT" 13 " , 

3303 IF (OXUM (KROW, KCOI ) +OXUM ( KROW , KCO) + .001) < (FUM(KRO.n, 

3900 IPRINT XPI ,";YPI.;" D" 

3 9 1 9 NEXT 

392 0 I PRINT” U " 

39 3'/ I PRINT” I 9 H 

3949 NEXT 

3950 REM PICT OX OTHER GRIP 
39 fir FOR KCOI~l TC 13 

397/ XI PO*XPO + 8i1* (KCOI -1 ) 

39G YI.PO-YPO + 15* (KCOI.-l ) 

3999 FOR KROW*l TO 29 
A 99/ KRO-KROW-1 
4919 IF KRO< 1 THEN KRO-1 
4929 XOFF*XI.PO-3H* (KROW-1 ) 

4939 YOFF-YI.PO+53* (KROW-1) 

4 (M 0 XPI * I NT (XOFF) 

A 9 5/ YPI * I NT (OXUM (KROW, KCOI ) * VSC AI E + YOFF+ • 5) 


KCOI ) +FUM (KROW, KCO) ) 


A 9 6 0 

I PRINT" 

1. 0 ” 

4 079 

IF (OXUM (KROW, KCOI ) +OXUM (KRO, KCOI ) + 

4 s o 

[.PRINT 

XPI;" , ” ; Y PI P ” 

4 9 99 

NEXT 


4100 

I. PRINT” 

U ” 

4110 

I PRINT” 

10 " 

A 1 2 0 

NEXT 


4 13 0 

[.PRINT” 

PH 290,0 n ” 

414 0 

RErt NO 

PI OT 

4 1 50 

END 



991 ) < ( rUM ( KROW , KCOI ) + FUM ( KRO , KC Of ) ) 


THEN 


THEN 


I n RI 


I n R I 


Figure B-2. (Concluded) 


R1/RD83-1 70 

B-12 


GAS-LIQUID DATA REDUCTION 


The gas-liquid data reduction process is somewhat more complex than the techniques 
used for liquid-liquid data reduction. The fluid simulants cannot be collected 
and measured directly in a sample grid as they are in the liquid-liquid collec- 
tion system since one component is a gas. The test procedure requires a probe 
that physically collects the liquid simulant (water), and senses the gaseous flow 
by total pressure, and gas sample measurements. A single gas liquid element flow- 
ing in a confined chamber induces very significant recirculation flow, which 
requires additional steps in the data acquisition and reduction. 

The use of a so called "base bleed" flow from the injection face area around the 
element, is an aid in supressing the recirculation flow. However, this gas flow 
becomes a portion of the gas flow in the element sampling area, and must be 
identified and separated from the computed flow. In our test setup, we accomplish 
the identification of the "base bleed" flow flux by using gaseous nitrogen as the 
element flow, and air as the "base bleed" flow. The probe gas flow is then 
sampled for oxygen content. The quantity of oxygen indicates the amount of bleed 
flow entrained in the sample. The data indicates a surprising amount of entrained 
flow, often being as much or more than the primary element flow (even in high mass 
flux zones near the center of the element zone). 

The test fixture allows probe movement along radial lines through the center of 
the sample chamber (a true rectangular coordinate system is not within the 
capabilities of the existing equipment), and the sample zone treatment in the data 
reduction program is based upon this radial system. The effective zone area is an 
annular arc section based on the probe radial distance from the center of the 
chamber, and the number of data "rays” indicated in the input data table. 
Figure B-3 shows this relationship for a case with eight rays assumed in the data 
acquisition. Frequently, in gas liquid testing, a limited number of actual sample 
"rays” are measured, and symmetry is assumed based on element type (anticipated 
shape of mass distribution). An extreme example of this technique is shown in 
Figure B-4. In this example, three actual "rays" have been measured, but an 
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equivalent of 12 rays are being used in the data reduci ion. The lense shaped mass 
distribution of the triplet dictates the type of symmetry used to establish the 
imaginary rays. 

The gas-liquid mixing data reduction program requires initial information on such 
items as the element configuration, flow rates, sampling time, radial spacing of 
the samples, and number of equivalent sampling rays. The individual probe sample 
point data is presented for each position, inputting the reference ray angle, the 
radial distance from the centerline, the dynamic pressure head reading, the amount 
of water collected during the sample time, and the percent of oxygen in the gases 
in the sample zone. 

From this data, the computer program calculates the value of oxidizer and fuel 
simulant flow flux in the sample area. The computation of liquid mass flow is 
very straight forward. The collected liquid mass divided by the sample time, is 
further divided by the probe open area to derive the liquid mass flux. The gas 
flux computation is less simple. The total pressure reading of the probe is 
related to gas velocity, but corrections must be made for the impact pressure con- 
tribution of the liquid droplets in the stream. The liquid flowrate is basically 
known, as is the gas density. The program uses a closed form computation based on 
the apparent fluid density of a two-phase flow of the gas and liquid. A further 
correction is made for the amount of gas from the base bleed. 

The two-phase flow correction, and the correction for entrained flow are rather 
large, and a potential source of error in this procedure. The entrained flow 
includes some recirculated element gases, as well as the "base bleed" flow, based 
on depressed oxygen content measured in zones well outside the element spray 
pattern." To correct for this effect, a depressed oxygen percentage (lower than 
the nominal air value) is input based on this outer zone measurement. This 
becomes the reference value for corrective gas flow flux values in the spray. 

One of the self-testing features of the data reduction procedure is the use of 
"collection efficiency." The program sums up the computed mass flux values times 
the respective sample areas, and compares this to the computed input flow values 
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for both the liquid and the gas. In many instances, collection efficiencies 
deviated significantly (by as much as 40%) from perfect collection efficiency. 
Since the probe measured areas do not overlap the entire sample zone (because of 
practical limits in the number of samples to be taken), some deviations between 
the computed total "collected” and the input values would be expected. However, 
the actual deviations in practice have been larger than anticipated, and several 
reviews of both the data acquisition procedures and data reduction procedures have 
been unable to resolve this. Since the liquid and gas flowrates are measured 
directly in these tests (using calibrated sonic orifices on the gas and a turbine 
meter on the liquid), it appears that some error in the test method exists. This 
problem still is being reviewed. Pending resolution of this problem, the error is 
assumed to be uniform across the entire sample area and the data is being 
"corrected” to the known flowrates. 

The data reduction program uses the corrected sample values to compute the classic 
mixing efficiency (E sub M) in the same way as the liquid-liquid mixing data pro- 
gram, and also computes the injection parameters, such as Elverum-Morey , momentum 
ratio, velocity ratio, velocity head ratio and the penetration parameter. The 
local sample injected mass flux values also are recorded in a data file for a 
graphic display on the plotter. In the plotted data shown in this report, only 
the actual "rays" where data was taken are plotted. The imaginary rays used to 
account for the total flowfield are not displayed. 

A listing of the gas/liquid data reduction program is shown in Fig. B-5. 
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1000 REM UPDATED VESION "CMXGLT" COI.D MIX GAS I.IQ TRIP 

1010 REM APPEND DATA STARTING AT I INE 9000 WITH MERGE COMMAND 

1020 REM UPDATED MAR 8 83 

1030 REM INCLUDES CORR FOR OX READING WITH BASE BLEED 
1,14 0 REM WITH OX CONTENT ADJUSTED FOR RECYCLED GN2 

1250 REM DATA REDUCTION PROGRAM FOR COLD FLOW LIQ MIXING VERSION WITH 
1060 REM CORRECTION FACTOR FOR MIXED FI OW DENSITY AT PROBE 
1073 PRINT" INPUT DAY AND TIME - ANY FORM" 

1080 INPUT DAYS 

1085 PRINT"TO SUPRESS AIL BUT SUMMARY PRINTOUT INPUT Y" 

1390 INPUT TPS 

1100 L PRINT" PROGRAM CMXGLT COLD FIOW MIXING DATA RED " ; DAYS 
1110 LPRINT 

1120 DIM WI.QS (40) ,WI QM (40) ,TS (4 0) f PC (4 0) ,WGS (4 0) , PRP (40) 

1133 DIM WGSF (40) ,WLQF(40) ,NZ(40) ,ZONA(40) ,WI.QZ(40) ,WGSZ(40) ,TEST(12) 

1 140 DIM OXPCTZ (40) 

1150 DIM RAY (40), RAD (40) 

1150 DIM CSTR (21) 

1170 DIM RGSFLX (40) ,RLQFLX (43) 

1183 REM 


1190 

REM C 

STAR 

DATA 




] 200 

DATA 

3000, 

3140, 

3290, 

3440, 

360 0 

1210 

DATA 

3950 , 

4130, 

4300, 

4470, 

4650 

1 220 

DATA 

4 900, 

5200, 

5700, 

6170, 

6200 

1230 

DATA 

5950, 

5400, 

4 500, 

3550, 

1 ooo 

124 3 

for r 

o 

H 

21 




1250 

READ 

CSTR (I ) 




1 2 fin 

NEXT 






200 

REM INPUT i 

CONS TAN 

TN 



2 010 

CAI IQ 

*1 




. 


2023 GNDNA*. 07274 
2033 GFUDN- 13.86 
2340 PDIA-.125 
2050 PAMB*1 3 . C 
2063 FLCD-.65 
2073 I. QDN *62.4 
2383 LOXDN»71 

2393 I.QSC-. 0022046 : REM CONV MI TO LBS 

2103 REM INITIAI IZE SUMMATION CONSTANTS 

2113 WL3CT-3 

2123 WGSCT-0 

2130 WTIQM-fl 

2143 WTGSM--C 

2150 EFMIX-1 

2163/ TEMFRC-3 

2170 ARSUM-0 

2180 DNGSM-0 

2190 CSMTOT-0 

2203 WTOTSM-O 

2210 REM 

3003 REM READ DATA 
3010 READ TESTS 

3020 READ DGSINJ,GSN,DLQINJ,I QN 

3333 REM DIA G5 INJ ORIF, NUMBER, DIA LQ INJ OR1F, NUMBER 
3340 READ KGSIN,WIQIN 


Figure B-5. Gas/Liquid Data Reduction Program 
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Ortt o 

OF PCG.i QUALITY 


3050 REM INJECTED FLOWRATES LB/SEC 
3060 READ NS,NRAY 

3070 REM NO DATA POINTS, NO OF RAYS TESTED 

3080 READ INPROB s REM INCREMENT BETWEEN PROBE RADIAL POS • N 

3090 READ TTIME.PCHB: REM TEST TIME AND CHBR PR 

3100 READ OXPB : REM OX PCT IN ENTRAINED CHAMBER GASES 

3110 REM SUGGEST HIGHEST OX VALUE MEASURED AT EDGE OF FAN 

3120 L PR I NT 

3130 LPRINT 

3140 LPRINT TESTS 

3150 LPRINT 

3160 LPRINT" PERCENT OX ASSUMED IN RECIRC GAS " ; OXPB 
3165 I.PRINT"TEST CHAMBER PRESSURE M ;PCHB 
3170 T PRINT 

3180 LPRINT" INPUT DATA AND INITIAL COMPUTATIONS" 

3190 LPRINT 

3192 LPRINT" INJ OR1F DIA - GAS " ; DCS I N J ; G 3N ; " PLCS - LIQ " ; DL QI N J ; L QN ; " PI 
3194 LPRINT" INJECTED WT FLOWS " ;WGSIN; " I B./SEC GAS - " ; WLCIN; “ I B/SEC I 10" 
3200 FOR N»1 TO NS 
3210 REM READ LOOP 

3220 READ NZ (N) : REM ZONE NUMBER FOR REF 

3230 READ RAY (N) , RAD (N) , PRP (N) , WLQS (N) , OX PCTZ (N ) 

3240 REM RAY ANGLE, RAD SAMP I.OC, PROBE PRESS , LIQ COt I Ml, PCT OX 
3250 IF TP$»"Y“ THEN 3280 

3260 I PRINT" ZONE" ; NZ (N) , " RAY" ; RAY (N) , ” RAD" ; RAD (N) /’PROBE PR" ; PRP (N 1 , 

3270 L PR I NT" LIQ ML" ;WLQS (N) , "PCT OX ” ; OX PCTZ ( N) 

3280 REM BYPASS PRINT 

3290 IF RAD(N)-0 THEN READ RAYMUI 

3300 IF PRP(N)-0 THEN PRP (N) ». 0000 0 1 

3310 ZONA (N) -RAD (N) * . 1 *2* 3 . 1 41 7/NRAY 

3320 IF RAD (N) < . 0001 THEN ZONA (N) » . 007854/NKAY 

3330 ZONA (N) -ZONA (N) *( INPROB/. 1 ) 

3340 ZONA (N) -ZONA(N) *RAYMUL 
3350 Wr QM (N) -WLQIN 
3 3 6u WGS ( N ) *WGS IN' 

3 3 7(4 TS(N)«TTIME 
3 380 PC (N ) * PCHB 

3390 WI QS (N) -WLQS (N) *LQSC/TS (N) : REM CORK SAMPLE FIOW 

340O IF WLQS (N) <.000001 THEN WLQS (N )-. 000001 
3410 ARSUM-ARSUM+ZONA ( N ) 

3420 GAMPN-NS 

3430 WTGSM-WTGSM+ (WGS (N) /SAMPN) 

3 4 4 0 WTI. QM - WTI. QM + ( wr QM ( N ) / S AM PN ) 

3450 NEXT N 

4000 REM 

4010 REM INITIAI MASS CALC I OOP 

4020 FOR N-I TO NS 

4030 REM GAS MASS Ft UX 

4040 QPSF-PRP (N) *144 

4050 GNDN-GNDNA* ( (PC (N) +PAMB) /I 4 . 7) 

4060 DNGSM-DNGSM+ (CNDN/SAMPN) 

4 070 WGSF(N)- ( (PRP(N) *GNDN/2.2 36) * .5) 

4080 PRINT"Q " ; QPSF ; " PSF FIRST WGAS " ; WGSF ( N ) ; ( WGSF (N ) *14 4) • " it/SEC ' GO 

4090 WGSF (N) -WGSF (N) * (WTGSM/WGS (N) ) 

4100 REM LIQ MASS FLUX 


Figure B-5. (Continued) 
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133 


4110 PAREA-.7854* (PDIA*2) 

4120 WLQF (N) *WI.QS (N) /PAREA 
4130 WLQFSF-WLQF(N) *144 

4140 PRINT" WT I.IQ FIOW PER SO FT/SEC" ;WIQFSF 
4150 WI.QZ (N)-WLQF(N) *20NA(N) 

4160 WI.QCT«WLQCT+WLQZ(N) 

4170 REM CORRECTION FOR MIXED FLOW DENSITY AT PROBE 
4180 CFSGS»WGSF(N)/GNDN 
4190 CFSCOL-WLQF(N)/60 
4200 CFST«CFSGS+CFSCOI. 

4210 WFST-WGSF (N) +WI.QF (N) 

4220 DNCOR-WFST/CFST 

4230 VHDMFI-(2.236*(WFST*2)/DNCOR) 

4240 VIIDFR-PRP (N) /VHDMFL 
4250 VHDEL-l-VHDFR 
4260 VHDCR-VHDEL/3 
4270 CORVHD-PRP(N) * (1-VIIDCR) 

4280 PRI NT" CORRECTED VEI HEAD (CORVHD*l 44) 

4 290 WGSF(N)- { ( (CORVIID*GNDN ) /2 . 2 36 ) “ . 5 ) 

4300 PR INT" NEK GAS FI 01* ";WGSF(N) ; (WGSF(N) *144) ; " ft/SEC/SQ FT" 

4310 KGSF(N)»WGSF (N) * ( WTGSM/WGS ( N ) ) 

4320 WGSZ (N) -WGSF (N) *ZONA (N) 

4330 REM NEW CORRECTION FOR TWO PHASE FLOW 
4340 A-1/(2*32.2*CNDN) 

4 350 B» (WI QFSF/64 .4) * ( ( 1/I.QDN ) + ( 1/GNDN ) ) 

4 36 0 C» ( (WI.?F5F“2) / (64 . 4*1 ODN) ) -QPSF 
4 370 X 1- (— 1 *B) - ( ( B* 2 ) -4 *A*C) “.5 
4080 X 1 -X 1 / ( 2 * A ) 

43.90 X2» (- 1 *B) + ( (B'2) -4*A*C) “ . 5 
4400 X2-X2/(2*A) 

44 10 PRINT" TWO SOI N ";X1;X2 

4420 PR I NT"GN' DENSITY " ; GNDN ; " EQUIV VEI ”;(X2/GNDN) 

4430 REM REPLACE GAS FLOW WITH NEW CORR 

4449 WGSF (N) -X2/144 

4450 REM CORR FOR HI EED AIR 

446 0 WGSF (N) -WGSF ( N ) * ( 1- (OXPCTZ ( N ) /CXPB ) ) 

4470 WGSZ (N) -WGSF (N) *70103 (N) 

4 4 00 PHI i GAS FIUX »" ;WGSF(N) GAS GOI I PER V. ONE " ;WGSZ(N) 

44: / REM CORR FOR NEC GAS FLUX 

44'. 4 IF WGSF (M) <. 00000] THEN WGSF ( N ) » . 0 0 0 0 3 1 
4490 W G SC T* WO SC T + WG S Z (N) 

4500 IF TP$»“Y" THEN 453C 

451!) I. PRINT" " , " G A S FIUX " ; W’GSF (N ) , " I IQ FIUX " ;WIQF (N) 

4520 PRINT 

4530 REM BYPASS PRINT 
4540 NEXT N 

4550 REM COLLECTION EFF 
4560 COI.EFC-WI QCT/WTLOM 
4570 COI EFO-WGSCT/WTGSM 
4580 I PRINT 

4590 l PRINT "COIIECTION EFF. I IQ*" ;COI EFC, "COLLECTION EFF. GS SIM«";COIEFO 

4600 I. PRINT 

4610 REM OVERALL O/F 

4620 RATIO-WTLQM/KTGSM 

4630 [.PRINT "OVERAII MIXTURE RATIO ( O/F) RATIO 


Figure B-5. (Continued) 
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4640 MEQCOR- ( (LOXDN/GFUDN) ' . 5) / ( (I.QDN/GNDN) * . 5) : HEM MOM RATIO CORR FOR EQ PROP 
4650 EQRAT-RATIO*MEQCOR 

4660 LPRINT "MOMENTUM EQUIVALENT MIXTURE RATIO (INJECTED) " ; EQRAT 
4670 LPRINT 

4680 LPRINT "MEAN FLOW RATES: " , "GO SIM-" ; WTGSM , " I. IQ-" ; WTLQM 

4690 LPRINT 

4700 WTOT-WTGSM+WTI QM 

4710 FRACT-WTGSM/'/TOT 

4 72 0 FI.XCM-WTLQM/ARSUM 

4730 FLXGSM-WTGSM/ARSUM 

4740 REM 

4750 IF TP$-"Y" THEN 4800 
4760 LPRINT 

4770 LPRINT " RAY" ," RAD" , “AREA" , "WT FLOW" , "MIX RAT.","EQIV 0/F","C STAR","G5 F!UX",“ 
4780 LPRINT"ANGI E" , " IN.","SQ IN" , "TOTAI " , " ","PROP." f " FPS" , " FRACT . “ , " FRACT" 

4790 LPRINT 

4800 REM BYPASS PRINT 

5000 REM FINAL SUMMATION LOOP 

5010 FOR N-l TO NS 

5020 REM CORRECT FOR COLI EFF 

5330 WLQZ (N) -WLQZ (N) /COI EFC 

5040 WGSZ (N) -WGSZ (N) /COI EFO 

5350 WTOTI -WLQZ (N) +WGSZ (N) 

5060 FRACTL-WGSZ (N) /WTOTL 

5070 IF FRACT! O FRACT THEN GOTO 5100 

5080 DECR- (WTOTL/WTOT) *{ ( FRACT-FRACTI. ) /( FRACT-1 ) ) 

5090 GOTO 5120 
5100 REM CONTINUE 

5110 DECR- (WTOTL/WTOT) * ( ( FRACT-FRACTI. ) / ( FRACT) ) 

5120 REM CONTINUE 

6133 EFMIX-EFMIX-DECR 

5143 REM NORMALIZED IOCAI MASS FLUX 

5150 RLQFI.X (N)-WLQZ (N)/(ZONA(N) *FI XCM) 

516 0 RGSFI.X (N) -WGSZ (N) / (ZONA(N) *FF XGSM) 

5170 RATL-WLQZ (N) /WGSZ (N) 

5180 EQRTI.-RATI "MEQCOR 
5190 EQOFRC-EQRTL/ ( EQRTL+ 1 ) 

5200 NRAT-INT (20*EQOFRC) 

5210 DRAT- ( 2 0*EQOFRC) -NR AT 
5220 CS1-CSTR (NRAT) 

5230 CS2-CSTR (NRAT+1 ) 

5243 DELCS-CS2-CS1 

5250 CSANS-CS1 +DRAT*DELCS 

5260 IF TP$-"Y" THEN 5280 

5270 IPRINT RAY (N ) , RAD (N ) ,ZONA(N) , WTOTL , RATI. , EQRTI , CSANS , RGSF! X (N) ,RIQFIX(N) 

5280 REM BYPASS PRINT 

5290 CSMTOT»CSMTOT+CSANS*KTOTI. 

5300 WTOTSM-WTOTSM+WTOTI 
5310 NEXT N 
5320 LPRINT 

5322 LPRINT"COI LECTED MIXTURE RATIO " ;WI.QCT/WGSCT 

5324 LPRINT" COMPUTED INJECTED 0/F " ; WTLQM/WTGSM 

5330 I. PR I NT "MOMENTUM EQI VALENT MIXTURE RATIO ( INJECTED) EQRAT 

5343 LPRINT "MIXING EFF.«";EFMIX 

5350 CSTRM-CSMTOT/WTOTSM 
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5360 LPRINT"MIXING LIMITED C STAR " ;CSTRM; " FPS" 

5370 REM O/AI-I. CSTAR 
5380 EQFROEQRAT/ (EQRAT + I ) 

5390 N RAT* I NT (20*EQFRC) 

5400 DRAT* (20*EQFRC ) -NRAT 
5410 CS1*CSTR (NRAT) 

5420 CS2-CSTR (NRAT+1 ) 

54 30 DEI.CS-CS2-CS1 

5440 CSTHE0-CS1+ (DRAT*DEI CS) 

54 50 I.PRINT M TI1E0 C STAR ";CSTHEC;" FPS" 

54 6 C NCSTR*1 00* (CSTRM/CGTIIEO) 

5470 I. PRINT "MIXING LIMITED C STAR EFF ".-NCGTR 
54 [PRINT 
54 0;'< REM 

rem injector operating conditions 
6010 GSAREA*GSN* (DGSINJ“2) *.7854 
60 2 0 GSVEI.-VvTGSM* 1 44 / ( DNGSM*GSAREA) 

6030 CFCOI *WTI QM/6:' 

6040 WTI Q-WTLQM 
6058 DNLQ-IQDN 

6060 LQAREA-I.QN* (DI.QINJ‘2) *.7854 
6870 I QVEL-WTI 0*144/ (DNI Q*I QAKEA) 

6080 RATMOM- (I QVEI *WTI.Q) / (G5VEL*VwTGSM) 

6398 AX* (WTGSM/WTLQ) “2 
6100 BX* (DNI Q/DNGSM) 

G1I0 CX- (I-QAREA/GSAREA) * 1 . 75 
6120 EI.MOR-AX*BX*CX 

6130 VrlDGS-3 . 62 5* (WTGSM~2) / (DNGSM* (GSN“2 ) * ( DGSINJ‘4 ) ■ 

614 0 VHDLQ-3.625* (WTI.Q'2) / (DNI.Q* ( LQN “ 3 ) * ( DI Q IN J * 4 ) ) 

6150 VHDR-VHDI.Q/VIIDGS 

6152 REM PENETRATION PARAMETER FOR GAS-IIQ-GAS TRIP 

6154 DIAOUT-DGSINJ 

6156 DI AIN»DI QINJ 

6158 DNOUT-DNCSM 

6160 DNIN-DNIQ 

6162 VEIOUT-GSVEl 

6164 VELIN-LQVEI 

6166 PEN 1-2. 5* (DI AOUT/DI A IN ) *( ( DNOUT/DN IN ) “ .5) 

6170 PEN2-.5* (VEI CUT/VEI IN) 

6172 REM .5 FOR 30 DEGREE ANGI E ( S IN 30) 

6180 PEN-PEN1 *PEN2 

6190 IPRINT "TRIPLET JET PENETRATION FACTOR*";PEN 
6200 L PR I NT 

6210 I. PR I NT "MOMENTUM RATIO-" ; RATMOM , " EL VERUM-MOREY COEFF-" ; El MOR 
6220 I.PRINT "INJECTION VEI- - GS SIM-" ;GSVEI , " I I QUID-" ; I QVEI 

6230 LPRINT "VEI HEAD - G S * " ; VHDGS ; " PSID - I. Q-" ; VHDI Q; " P5ID -VEI HEAD RATIO*";VHL 

6240 LPRINT 

6250 LPRINT 

6260 LPRINT 

6270 LPRINT 

7000 REM 

7010 REM FILE FOR PLOT 

7020 PRINT" INPUT Y IF PLOT FILE IS DESIRED" 

7030 INPUT FICOD5 

7040 IF FICODSO "Y" THEN 8999 
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7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7150 

7160 

7170 

7180 

7190 

7200 

7205 

7210 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7420 

74 30 

744 P 

7450 

74 6 0 

0999 

9000 

9010 

9020 

9030 

9035 

9040 

9050 

9100 

9110 

9120 

9130 

9140 

9150 

9160 


QRIQ.K; X 

OF FOGR 




wi 


, ;CHR${34) ; TESTS ;CHR$ (34) 
” ; EQRAT 
" ; EFM I X 


HAY 


RAY (NSTRT) 


PRINT" INPUT DESIRED NAME OF FILE" 

INPUT NFI $ 

NFI$*NFI $+" .BAS" 

OPEN" 0" , # 1 , NF I $ 

REM 

PRINT# 1 ,"10 DATA 
PRINT# 1 , " 20 DATA 
PRINT # 1 , " 3 0 DATA 
NSTRT-1 
FOR L»2 TO NS 

IF L*NS THEN L»L+1:G0T0 7180 
IF RAD (L)>.0001 THEN 7400 
REM LAST SET OF DATA 
P K I NT # 1 , I * 1 0 0 ; " REM DA TA FOR 
PRINT# 1 , (1*100 + 1) ; "DATA " ; 

NDT*1 + ( L -NSTRT) 

WRITE# 1 , RAY ( NSTRT) , NDT 

FOR N-NSTRT TO L 

PR I NT # 1 , (T *100+ ( 2 * N ) ) ; "DATA 

IF N»L THEN 7280 

WRITE# 1 , RAD (N) , R[ QFI X (N) 

GOTO 7290 

W R I T E # 1 , (HAD(N-1)+. 1) ,0 

REM 

NEXT 

FOR N* NSTRT TO L 

PR I NT # 1 , (1*100 + 50+ (2*N) ) ; "DATA " ; 

IF N-L THEN 7360 
WRITE# 1 , RAD (N) , RGSFLX ( N ) 

GOTO 7370 

WRITE » 1 , (RAD (N-l ) + . 1 ) , C 

REM 

NEXT 

NSTRT*! 

REM 
REM 
NEXT 
CT OSE ( 1 ) 

END 

REM 

REM LAST IINE BEFORE DATA 

DATA "TRIPIKT NO 2, NOM MOM R, 50 PSIG, I IN" 

• 063 f 2, . 04 5, 1:REM GAS INJ ORIF DIA AND NO. I IQ ETC” 
.0138, .042: REM GAS WT FI OW I IQ a T FLOW 

12, 8: REM NO OF DATA POINTS, NO OF RAYS 

. 1 : R EM INCREMENT BETWEEN SAMPLES 
60 , 50 : REM TEST TIME IN SEC, TANK PRESS PSIG 
16.5 : REM PERCENT OX TO BE USED TO FACTOR OUT 

2.85, 33, 10.5 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


RECTRC 


0, 0, 

6 

2 


■ 1» 1.45, 30, 11 
>2, .25, 8, 12.5 
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9170 

DATA 

4 




9180 

DATA 

0, 

.3, 

0. 1, 

14.5 

920 0 

REM 





9210 

DATA 

5 




9220 

DATA 

90, 

0, 

2.7, 

40, 11 

9230 

DATA 

2 




9240 

DATA 

6 




9250 

DATA 

90, 

.1, 

2.55 

, 40, 11 

9260 

DATA 

7 




9280 

DATA 

90, 

• 2, 

1.7, 

36, 11.5 

9290 

DATA 

8 




9300 

DATA 

90, 

.3, 

.85, 

32, 13 

9310 

DATA 

9 




9320 

DATA 

90, 

.4, 

.4, 

31, 13.5 

9340 

DATA 

10 




9350 

DATA 

90, 

• 5, 

• 25, 

26, 14.5 

9360 

DATA 

11 




9370 

DATA 

90, 

.6, 

.1, 

21 , 15 

9380 

DATA 

12 




9390 

DATA 

90, 

-7, 

.05, 

8, 16 

94 0 0 

DATA 

13 




94 10 

DATA 

90 , 

.8, 

0, 5 

, 16 

9999 

END 
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